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INTRODUCTION. 

§ 1. A great many if not all of the problems in mathematics 
may be so formulated that they consist in finding from given 
data the values of certain unknown quantities subject to certain 
conditions. We may distinguish different stages in the solution 
of a problem. The first stage we might say is the proof that the 
quantities sought for really exist, that it is possible to satisfy 
the given conditions or, as the case may be, the proof that it is 
impossible. In the latter case we have done with the problem. 
Take for instance the celebrated question of the squaring of the 
circle. We may in a more generalized form state it thus: Find 
the integral numbers, which are the coefficients of an algebraic 
equation, of which t is one of the roots. Thirty years ago 
Lindemann showed that integral numbers subject to these con- 
ditions do not exist and thus a problem as old almost as 
human history came to an end. Or to give another instance 
take Fermat's problem, for the solution of which the late Mr. 
Wolfskehl, of Darmstadt, has left $25,000 in his will. Find the 
integral numbers x, y, z that satisfy the equation 

ar" + y* = 2*, 

where n is an integral number greater than two. Fermat main- 
tained that it is impossible to satisfy these conditions and he is 
probably right. But as yet it has not been shown. So the 
solution of the problem may or may not end in its first stage. 
In many other cases the first stage of the solution may be so 
easy, that we immediately pass on to the second stage of finding 
methods to calculate the unknown quantities sought for. Or 
even if the first stage of the solution is not so easy, it may be 
expedient to pass on to the second stage. For if we succeed in 
finding methods of calculation that determine the unknown quan- 
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VI GRAPHICAL METHODS. 

tities, the proof of their eidstence is included. If on the other 
hand, we do not succeed, then it will be time enough to return 
to the first stage. 

There are not a small number of men who believe the task of 
the mathematician to end here. This, I think, is due to the 
fact that the pure mathematician as a rule is not in the habit of 
pushing his investigaticfti so far as to find something out about the 
real things of this worid. He leaves that to the astronomer, to 
the physicist, to the engineer. These men, on the other hand, 
take the greatest interest in the actual numerical values that 
are the outcome of the mathematical methods of calculation. ' 
They have to carry out the calculation and as soon as they do so, 
the question arises whether they could not get at the same result 
in a shorter way, with less trouble. Suppose the mathematician 
gives them a method of calculation perfectly logical and con- 
clusive but taking 200 years of incessant numerical work to 
complete. They woidd be justified in thinking that this is not 
much better than no method at all. So there arises a third stage 
of the solution of a mathematical problem in which the object is 
to develop methods for finding the result with as little trouble as 
possible. I maintain that this third stage is just as much a 
chapter of mathematics as the first two stages and it will not do 
to leave it to the astronomer, to the physicist, to the engineer or 
whoever applies mathematical methods, for this reason that 
these men are bent on the results and therefore they will be apt 
to overlook the full generality of the methods they happen to 
hit on, while in the hands of the mathematician the methods 
would be developed from a higher standpoint and their bearing 
on other problems in other scientific inquiries would be more 
likely to receive the proper attention. 

The state of affairs today is such that in a number of cases the 
methods of the engineer or the surveyor are not known to the 
astronomer or the physicist, or vice versa, although their prob- 
lems may be mathematically almost identical. It is particularly 
so with graphical methods, that have been invented for definite 
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problems. A more general exposition makes them applicable 
to a vast number of cases that were originally not thought of. 
In this course I shall review the graphical methods from a 
general standpoint, that is, I shall try to formulate and to teach 
them in their most generalized form so as to facilitate their 
application in any problem, with which they are mathematically 
connected.^ The student is advised to do practical exercises. 
Nothing but the repeated application of the methods will give 
him the whole grasp of the subject. For it is not sufficient to 
understand the underlying ideas, it is also necessary to acquire a 
certain facility in applying them. You might as well try to learn 
piano playing only by attending concerts as to learn the 
graphical methods only through lectures. 

* For the literatiire of the subject see '^Encyklopfidie der mathemalischen 
Wissenschafteii/' Art. R. Mehmke, '' Numerisches Rechnen," and Art. F 
WUleiB and C. Runge, "Graphische Int^;ration." 



CHAPTER I. 
Graphical Calculation. 

§ 2. Graphical Arithmetic. — ^Any quantity susceptible of mensu- 
ration can be graphically represented by a straight line, the 
length of the line corresponding to the value of the quantity. 
But this is by no means the only possible way. A quantity 
might also be and is sometimes graphically represented by an 
angle or by the length of a curved line or by the area of a square 
or triangle or any other figure or by the anharmonic ratio of four 
points in a straight line or in a variety of other ways. The 
representation by straight lines has some advantages over the 
others, mainly on account of the facility with which the ele- 
mentary mathematical operations can be carried out. 

What is the use of representing quantities on paper? It is a 
convenient way of placing them before our eye, of comparing 
them, of handling them. If pencil and paper were not as cheap 
as they are, or if to draw a line were a long and tedious under- 
taking, or if our eye were not as skillful and expert an assistant, 
graphical methods would lose much of their significance. Or, 
on the other hand, if electric currents or any other measurable 
quantities were as cheaply and conveniently produced in any 
desired degree and added, subtracted, multiplied and divided 
with equal facility, it might be profitable to use them for the 
representation of any other measiuable quantities, not so easily 
produced or handled. 

The addition of two positive quantities represented by straight 
lines of given length is effected by laying them off in the same 
direction, one behind the other. The direction gives each line a 
beginning and an end. The beginning of the second line has to 
coincide with the end of the first, and the resulting line represent- 
ing the sum of the two runs from the beginning of the first to 
2 1 
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the end of the second. Similarly the subtraction of one positive 
quantity from another is effected by giving the lines opposite direc- 
tions and letting the beginning of the line that is to be subtracted 
coincide with the end of the other. The result of the subtrac- 
tion is represented by the line that runs from the beginning of 
the minuend to the end of the subtrahend. The result is positive 
when this direction coincides with that of the minuend^ and nega- 
tive when it coincides with that of the subtrahend. This leads 
to the representation of positive and negative quantities by lines of 
opposite direction. The subtraction of one positive quantity from 
another may then be looked upon as the addition of a positive and 
a negative quantity. I do not want to dwell on the logical explana- 
tion of this subject, but I want to point out the practical method 
used for adding a large number of positive and negative quantities 
represented by straight lines of opposite direction. Take a 
straight edge, say a piece of paper folded over so as to form a 
straight edge, mark a point on it, and assign one of the two 
directions as the positive one. Lay the edge in succession over 
the different lines and run a pointer along it through an amount 
equal in each case to the length of the line and in the positive 
or negative direction according to the sign of the quantity. The 
pointer is to begin at the point marked. The line running from 
this point to where the pointer stops represents the sum of the 

given quantities. The advan- 
tage of this method is that the 
intermediate positions of the 
pointer need not be marked pro- 
vided only that the pointer keeps 
its position duritig the move- 
ment of the edge from one line 
to the next. As an example take 
the area. Fig. 1. A number of 
rectangular strips ^ cm. wide are substituted for the area so that, 
measured in square centimeters, it is equal to half the sum of 
the lengths of the strips measured in centimeters. The straight 




Fio 1. 



GBAPmCAL CALCULATION. 3 

edge is placed over the strips in succession and the pointer is 
run along them. The edge is supposed to cany a centime- 
ter scale and the pointer is to begin at zero. The final position 
of the pointer gives hfdf the value of the area in square centi- 
meters. The drawing of the strips may be dispensed with, their 
lengths being estimated, only their width must be shown. If 
the scale should be too short for the whole length, the only thing 
we have to do is to break any of the lengths that range over the 
end of the scale and to count how many times we have gone 
over the whole scale. I have found it convenient to use a little 
pointer of paper fastened on the runner of a slide rule so that it 
can be moved up and down the metrical scale on one side of the 
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slide rule. The area is in this manner determined rapidly and 
with considerable accuracy, very well comparable to the ac- 
curacy of a good planimeter. If the area of any closed curve 
is to be found, the way to proceed is to choose two parallel 
lines that cut off two segments on either side (see Fig. 2), to 
measure the area between them by the method described above 
and to estimate the two segments separately. If the curves of 
the segments may with sufficient acciuracy be regarded as arcs 
of parabolas the area would be two thirds the product of length 
and width. If not they would have to be estimated by substitute 
ing a rectangle or a number of rectangles for them. 
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In the same way the addition and subtraction of pure numbeis 
may also be carried out. We need only represent the numbers 
by the ratios of the lengths of straight lines to a certain fixed 
line. The ratio of the length of the sum of the lines to the length 
of the fixed lines is equal to the sum of the numbers. The con- 
struction also applies to positive and negative nimibers, if we 
represent them by the ratio of the length of straight lines of 
opposite directions to the length of a fixed line. 

In order to multiply a given quantity c by a given number, 
let the number be given as the ratio of the lengths of two straight 
lines a/b. If the quantity c is also represented by a straight line, 
all we have to do is to find a straight line x whose length is to 
the length of c as a to b. This can be done in many ways by 





Fig. 3. 



Fig. 4. 



constructing any triangle with two sides equal to a and 6 and 
drawing a similar triangle with the side that corresponds to 6 made 
equal to c. As a rule it is convenient to draw a and b at right 
angles and the similar triangle either with its hypotenuse parallel 
(Fig. 3) or at right angles (Fig. 4) to the hypotenuse of the first 
triangle. Division by a given number is effected by the same con- 
struction; for the multiplication by the ratio afb is equivalent 
to the divisions by the ratio 6/a. 

If a, 6, c are any given numbers, we can represent them by the 
ratios of three straight lines to a fixed line* Then the ratio of 
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the line constructed in the way shown in Fig. 3 and Fig. 4 to 
the fixed line is equal to the number 

ac 
6* 

MuItipUcation and division are in this way carried out simul- 
taneously. In order to have multiplication alone, we need only 
make h equal 1 and in order to have division alone, we need only 
make a or c equal 1. 

In order to include the multiplication and division of positive 
and negative numbers we can proceed in the following way. Let 
the lines corresponding to a, x, Fig. 3, be drawn to the right side 
of the vertex to signify positive numbers and to the left side to 
signify negative numbers. Similarly let the lines corresponding 
to by c be drawn upward to ^gnify positive numbers and down- 
ward to signify negative numbers. Then the drawing of a 
parallel to the hypotenuse of the rectangular triangle a, 6 through 
the end of the line corresponding to e will always lead to the 

number 

oc 

whatever the signs of a, 6, c may be. 

The same definition will not hold for the construction of Fig. 4. 
If the positive direction of the line corresponding to a is to the 
right and the positive direction of the line corresponding to b is 
upwards then the positive directions of x and c ought to be such 
that when the right-angled triangle x, c is turned through an 
angle of 90^ to make the positive direction of x coincident 
with the positive direction of a, the positive direction of c coin- 
cides with the positive direction of 6. If we wish to have the 
positive direction of x upward, the positive direction of c would 
have to be to the left, or if we wish to have the positive direction 
of e to the right, the positive direction of x would have to be 
downward. If this is adhered to, the construction for division 
and multiplication will include the signs. 
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§ 3. Integral Functions. — ^We have shown how to add, subtract, 
multiply, divide given numbers graphically by representing them 
as ratios of the lengths of straight lines to the length of a fixed 
line and finding the result of the operation as the ratio of the 
length of a certain line to the same fixed line. By repeating 
these constructions we are now enabled to find the value of any 
algebraical expression built up by these four operations in any 
succession and repetition. Let us see for instance how the values 
of an integral function of x, that is to say, an expression of the form 



may be 



found by geometrical construction, where ao, ai "' anyX 

are any positive or negative 
numbers. We shall first as- 
sume that all the numbers are 
positive, but there is not the 
least difficulty in extending 
the method to the more gen- 
eral case. 

Now let Oo, ai, fh, '" an 
signify straight lines laid off 
on a vertical line that we call 
the y-axis, one after the other 
as if to find the straight line 
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00+01 + 02+ ••• + o«. 



The lengths of these lines measured in a conveniently chosen 
unit of length are equal to the numbers designated by the same 
letters. In Pig. 5 oo runs from the point to point Ci, oi from 
CitoCt, '"an from Cn to C«+i. 

Let a; be the ratio of the lines Ox and 01, Fig. 5, drawn hori- 
zontally from to the right. The length 01 is chosen of con- 
venient size independent of the unit of length that measures the 
lines Oo, Oi, • • • an. The length Ox is then defined by the value 
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of the ratio x. Through x and 1 draw Unes parallel to the 2^-axis. 
Through Cn+i draw a line parallel to Ox, that intersects those two 
parallds in Pn and Bn- Draw the lineBn^n that intersects the 
parallel through x in Pn~i. Then the height of Pn~t above Cn 
will be equal to OnX. For if we draw a line through Pn^i parallel 
to Ox intersecting the j^-axis in D^ the triangle CnDnPn-i will be 
similar to CnCn^iBn and their ratio is equal to x, therefore 
CnDn = anX. Couscquently the height of P^-i above Cn^i is 
equal to Cn^iDn = anX + an^i. Now let us repeat the same 
operation in letting the point Pn take the part of Cn+i. Through 
Dn draw a line parallel to Ox, that intersects the parallels through 
X and 1 in Pn^i and Bn^i. Draw the line Bn^iCn~i that intersects 
the parallel through x in Pn^» 
Then the height of Pn^ above 
Cn-i will be equal to 

Cn^lDn'X = (a^ + an^x)x, 

and the height above (7»-s will be 
equal to 

Ana!* + cifi^iX + tti^-a. 

Continue in the same way. Draw 
Pnr^n^ parallel to Ox, draw 
Bn^Cn^ and find the point Pi»-s. 
Then the height of P»-i above C 
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will be 



{an^ + an^ix + an^)x 
and the height of P«^ above (7«-4 

Finally a point Po is found (see Fig. 6 f orn = 4) by the inter- 
section of BiCi with the parallel to the 2^axis through x, whose 
height above is equal to 

anX"" + an^ia:*"^ + l-oix + oo 

Let us defflgnate the line xPq by y, so that 
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in the sense that y b a vertical line of the same direction and 

length as the sum of the vertical lines a^^ a^ ig*^» * * * aio;, oo. 

The same construction hdds good for values of x greater than 

1 or nq^tive. The only difference is that the point x is beyond 

the interval 01 to the right of 1 or to the left of 0. The negative 

sign of 

a»a:, a^ + a«^i, a^ + a«_ix, etc, 

will signify that the direction of the lines is downward. Nor are 
any alterations necessary in order to include the case that several 
or all of the lines Oo, oi, • * • €Ln are directed downward and corre- 
spond to nq^tive numbers. They are laid off on the y-azis in 
the same way as if to find the sum 

Oo + ai + Oi + ••• '\' €Ln9 

CmW lying above or below (7« according to a^ being directed 
upward or downward. The construction can be repeated for a 
number of values of x. The points Po will then represent the 
curve, whose equation is 

y = Oo + aia: + h a»a:», 

X and y measuring abscissa and ordinates in independent units 
of length. 

In order to draw the curve for large values of a; a modification 
must be introduced. It will not do to choose 01 small in order 
to keep X on your drawing board; for then the lines B«0« will 
become too short and thus their direction will be badly defined. 
The way to proceed is to change the variable. Write for instance 
X ^ xflOf so that X is ten times as small as x and write 

-4. = a.- 10*. 
Then as 
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we find 

Lay off the lines Aq, Ai, • • • An in b, convenient scale and let 
X play the part that x played before. The curve differs in scale 
from the first curve and the reduction of scale may be different 
for abscissas and ordinates but may if we choose be made the 
same so that it is geometrically similar to the first curve reduced 
to one tenth. It is evident that any other reduction can be 
effected in the same manner. By increasing the ratio x/X we 
enhance the value of ^» in comparison to the coefficients of lesser 
index, so that for the figure of the curve drawn in a very small 
scale all the terms will be insignificant except AnX\ In this 
case the points Ci, C%, • • •, Cn will very nearly coincide with 
and only Cn+i will stand out. 

It is interesting to observe that the best way of calculating an 
integral function 

Oo + aiar + a^ + • • • + OnX* 

for any value of x proceeds on exactly the same lines as the 
geometrical construction. The coefficient an is first multiplied 
with X and an^i is added Call the result On^i'. This is again 
multiplied by x and an^s is added. Call this result an^\ Con- 
tinuing in this way we finally obtain a value of (io\ which is equal 
to the value of the integral function for the value of x considered. 
Using a slide rule all the multiplications with x can be effected 
with a single setting of the instrument. The coefficients a« and 
the values ctj are best written in rows in this way 

Gn ttn^i a«-8 • • • tti Oq 

anX an^iX ' ' ' djX a\X 

an^i an-Jk '" On! Oo' 

The accuracy of the slide rule is very nearly the same as the 
accuracy of a good drawing. But the rapidity is very much 
greater. When therefore only a few values of the integral func- 
tion are required, the geometrical construction will not repay 
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the trouble. It is different, however, when the object is to make 
a drawing of the curve. The values supplied by calculation 
would have to be plotted, while the geometrical construction 
furnishes the points of the curve right away and in this manner 
gains on the numerical method. 

There is another geometrical method, which in some cases 
may be just as good. Let us propose to find the value of an 
integral function of the fourth degree. 

y == Oo + aiir + OfiO? + a^ + a^ 

and let all coefficients in the first instance be positive. 

The coefficients Osh ^^ly as> os, 04 are supposed to be represented 
by straight lines, while x will be the ratio of two lines. The lines 

<>09 Aiy os, 08, 04 are laid off in a 
broken line ao to the right from 
Co to (7i, ai upward from C\ to 
C2, 02 to the left from C% to 0%^ Os 
downward from 0% to C4, 04 again 
to the right from C4 to Cj (Fig. 7). 
Through Cs draw a line 0^4 to 
a point A on C%C^ or its prolonga^ 
tion and let z be equal to the 
ratio C^ : C4C5 taken positive 
when C42I has the same direc- 
tion as CfCi. Then we have 




and 



C^ SB Otf?, 

CtA = otf? + oj. 



C^ and CtA are positive or negative according to their direction, 
being the same as the direction of C%C^ or opposite to it. Through 
A draw the line AB forming a right angle with C^ to a point B 
on CtC% or its prolongation. Then we have 



and 



CsjB = C%A*z = (04a; + oj) a; 
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CtB = 04Z" + 0|X + dj. 

CiB and CtB are positive or negative according to thdr direction 
being the same as the direction of CtCt or oppoate to it Kmi- 
larly we get 

CiD = 041* + otr* + djx + Oi. 
and finally 

CoB = a^a:* + a^ + a^ + oia: + oo. 

CqE is poative, when £ is on the right side of Co and negative 

when on the left side. When the point A jnoves along the line 

CtCi, the prant E will move 

along the line CaCi and its 

position will determine the 

values of the integral function. 

To find the portion of E for 

any position of ^, we might 

use transparent squared paper, 

that we pin onto the drawing 

at Ci, 80 that it can freely be 

turned round C^ Following 

the lines of the squared paper 

along ChABDE after turning it through a small an^e fumishea 

the position of £ for a new position of A (Fig. 8). 

To indiide the case of negative coefiSdents we draw the corre- 
sponding line in the opposite direction. If for instance 04 is 
negative CtC^ would have to lie above Ct; but CtA would have 
to be counted in the same way as before, po^tive in a downward, 
negative in an upward direction. 

The ezten^on of the method to integral functions of any degree 
is obvious and need not be inmsted on. It may be applied with 
advantage to find the real roots of an equation of any degree. 
For this purpcMie the broken line C^BDE would have to be 
drawn in such a way that E coinddes with Co. In the case of 
Tig. 7, for instance, it is easily seen that no real root exists. 
Hg. 9 sibows the application to the quadratic equation. A drde 
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C» 



is drawn over CoCi as diameter* Its intersections with CiCt 
furnish the points A and A' that correspond to the two roots. 
Both roots are negative in this case. 

The first method of constructing 
the values of an integral function can 
be extended to the case where the 
function is given as the sum of a 
number of polynomials of the form 

y = oo + ai(x —p)+ (h{x—p)(x —9) 
+ az{x — p){x — q){x — r) + • • •. 

Let us again suppose ao» a^ os^ ' * * 
Pjo^ g^ to represent straight lines laid off as 

before on the y-axis upwards or down- 
wards as if to find their sum. x, p^Qft • • - are meant, to be num- 
bers represented by the ratio of certain segments on the axis of 
abscissas. Let us consider the case of f oiu* terms, the highest poly- 
nomial being of the third degree. The fixed distance between the 
points marked p and p+ l,q and q+ 1, r and r + 1 on the 
axis of abscissas, Fig. 10 is chosen arbitrarily and the position 
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of the points marked p, q, r, x is made such that the ratio of 
Op, Oq, Or, Ox to that fixed distance is equal to the numbers 
p, q, r, X. For negative values the points are taken on the left 
ofO. 
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Draw parallels to the y^sjis through p, q, r, x, p+l, q+ 1, 
r+h On the parallel through r+1 find the point Qo of the 
same ordinate as C4 and on the parallel through r find the point 
Aq of the same ordinate as C|. Join Aq and Qo by a straight 
line and find its intersection Pi or that of its prolongation with 
the parallel through x. The height of Pi above Ci or Ao is 
equal to <h(x — r) and the height above C2 is equal to ozix — r) 
+ a%. On the parallel through 9 + 1 find a point Qi of the same 
ordinate as Pi and on the parallel through q a point Ai of the 
same ordinate as Cs. Join Ai and Qi by a straight line and find 
its intersection Ps or that of its prolongation with the parallel 
through X. The height of P2 above Cs or Ai is equal to 

[at(x - r) + (h](x - g), 

and the height above Ci is equal to 

(h{x — r)(x — q) + (h(x — q) + (h. 

Finally find a point Q2 on the parallel through p + 1 of the 
same ordinate as P2 and a point A2 on the parallel throiigh p of 
the same ordinate as Ci. Join A2 and Q2 by a straight line and 
find its intersection Pz or that of its prolongation with the par- 
allel through X. The height of P| above Ci or At will then be 

equal to 

[Oiix - r){x — q) + (h{x — q) + ai](x — p) 

and the ordinate of Pz will be eqiial to the given integral function 

y == ozix— r){x — q)(x — p) + (h{x — q){x — p) 

+ ai(x — p) + oo. 

For large numbers p, q, r, x we use a similar device as before by 
introducing new numbers P, Q, i2, X equal to one tenth, or one 
hundredth or any other fraction of pqrx. For instance 

P = p/10, (2 = ?/10, fl=r/10 X^mO. 

We then write 

Ao =00, Ai^ lOoi, At = 10002, Az = lOOOoi, 
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and obtain 

y ^ A. + Ai(X - P) + MX - P)a - (R 

+ ^(Z-P)(Z-Q)(Z-il). 

The scale for tbe fines A^, Ay, At, A% and y mnst then be reduced 
conveniently and the values are constmcted in the same way as 
before. 

Now let us consider the inverse problem. The values of the 
integral function are given for 

a? = P, ?, ff «; 

find the fines Oo, Oi, Oif Os, so that the value dt the int^r&I function 
may be found for any other value iAxm the way shown above. 

Let us designate the given values of the integral function for 
^^ Vf^f^f^ by Vv* y^f Vrf y* c^d the points on the parallels through 
p, q^ r, 9 with these ordinates by P,Q,R,S (see Fig. 12). 

For:c s p the integral function 

y- ao + (h(x — p) + (h(x — p)(ar — q) + a^ix— p){x—q)(x—r) 

reduces to oq. Therefore we have j^p = oo- The point Ci is 
found by drawing a parallel to the axis of abscissas through P 

and taking its intersection with 
the axis of ordinates. 

In order to find Ct draw a 
straight line through P and Q 
and find its intersection A with 
the parallel through p + 1 (Fig. 
11). A parallel to the axis of 
abscissas through A intersects 
the axis of ordinates in Ct. For 
the differences yq—yp and ya^yp 
(writing ya for the ordinate of 
A) are proportional to the differences of the abscissas and con- 
sequently in the ratio (g — p) : 1. Therefore 
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In the same way as the point Q on the parallel through q we 
might join any point Z on a parallel through x with the point P, 
find the intersection with the parallel through p + 1 and draw a 
parallel to the axis of absdssas. The point of intersection of 



this parallel with the vertical through x let us call X' and its 
ordinate j/. Then we have 

y'-yp = rir?= «» + °*(^ - ?) + ««(* - ?)(* - r). 

X p 

Let us carry out this construction not only tat x •= q but also 
for X = r and x = s. This leads us to three points Q', R', S' 
on the verticals through q, r, s, whose ordinates are the values 
of the integral functions 

y" = (oo + oi) + otix — q) + OjC* — q){x — r). 
In this way we have reduced our problem. Instead of having 
to find an integral function of the third degree from four given 
points P, Q, R, S, we have now only to find an integral function 
of the second degree from three given points Q', R', S'. A second 
reduction is effected in exactly the same manner. Q' is joined 
with R' and 8' by straight lines and through their intersection 
with the vertical through q+ 1 paraUels to the axis of abscissas 
are drawn that intersect the verticals through r and « in the 
points R" and S" respectively. The ordinates of these points 
are the values of the integral function y" defined by 
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lor z— r and a: = *, or 

y" = ao + ai + oj + (h{x — r). 

The horizontal through R" intersects the axis of ordinates in the 
point Cz. Finally we find C4 by drawing a parallel to the axis 
of abscissas through the intersection of R'S" or its prolongation 
with the vertical through r + 1. 

Having found the points CiC^CzCi we can now for any value 
of X construct the ordinate 

y = ao + ai{x — p) + 02(3: — p)(x — q) 

+ (h(x — pXx - q){x - r), 

and thus draw the parabola of the third degree passing through 
the four points P, Q, R, S. 

The construction may be somewhat simplified first by making 
p + 1 = g. Our data are the points P, Q, -R, S, and we are 
perfectly at liberty to make the vertical through p + 1 coincide 
with the vertical through Q. In this case the point Q' will 
coincide with Q. The parabola of the second degree through the 
points Q'iJ'S' is again independent of the distance between the 
verticals through q and 9 + 1 and at the same time independent 
of the point P. Therefore we are perfectly at liberty, for the 
construction of any point of this parabola, to make the vertical 
through q+ 1 coincide with the vertical through R even if the 
distance of the verticals through P and Q is different from that 
of the verticals through Q and -R. R" will in this case coincide 
with R\ The procedure is shown in Fig. 12. Starting from 
the points P, Q, R, S the first step is to find -R', S' by connecting 
R and 8 with P and drawing horizontals through the inter- 
sections Ar, Aa with the vertical through q. The next step is to 
find /S" by connecting Q (identical with Q') with S' and drawing 
a horizontal through the intersection with the vertical through r. 
Now the straight line JB"S" can be drawn (-R" being identical 
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with R'). On the vertical through any point x take the inter- 
section with R"S" and pass horizontally to the point Ax' on the 
vertical through r. Draw the line Q'Ax' and find its intersection 
with the vertical through x. This point is on the parabola 
through Q'R'S'. Pass horizontally to the point As on the 
vertical through q and draw the line AgP. Its intersection with 
the vertical through ;r is a point on the parabola of the third 
degree through P, Q, R, S. 

The method is evidently applicable to any number of given 
points, the degree of the parabola being one unit less than the 
number of points. 

The methods for the construction of the values of an integral 
function may be applied to find the value of any rational function 

y = R(x). 

For a rational function can always be reduced to the form of a 
quotient of two integral fimctions 

Rix) = gi(x)/g,{x). 

Now after having constructed curves whose ordinates give the 
values of gi(x) and gt{x) for any abscissa x (Fig. 13), R(x) is found 
in the following manner. 
Through a point P on the iirAzia 

axis of abscissas draw a 
parallel to the axis of or- 
dinates. Let Ox and 0% 
be the points whose ordi- 
nates are equal to gi{x) 
and g%(x). Pass horizon- 
tally from Oi to Oi on the 
vertical through P and 
from Qt to Ot on the axis of ordinates. Draw a line through 
P and 6t and produce it as far as ^ where it intersects the 
horizontal through ft. Then R(x) is equal to the ratio OiA 
to PO. Oi'A may then be set off as ordinate on the vertical 
3 
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through z and defines the point M whose ordinate is equal to 
R{x) in lengthy when OP is chosen as the unit of length. 

§4. Linear Functions of Any Number of Variables. — ^Let us 
consider a linear function of a number of variables xi, oci *•• Xn, 

Oo + a^Xi + a^Xi + • • • + an;Xn, 

where oo, oi^ 02, * • * an are given numbers positive or negative. 
The question is how the value of this linear function may be 
conveniently constructed for various systems Xi, Xi, ••• Xn. 
Suppose Oo, ai, • • • an to represent horizontal lines directed to 
the right or left according to the sign of the corresponding number 
and to be laid off on an horizontal axis in succession as if to find 
the sum 

00+^1+08+ •" + Onf 

ao begins at and nms to Cu W begins at C\ and runs to C2 and 
so on (Fig. 14). The numbers Xi^ Xg^ • • • 2:» let us represent 
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by ratios of lengths. We draw a vertical line through and 
choose a point P on the horizontal axis. Then let xi be equal 
to the ratio 01/PO, X2 = 02/PO, etc. If P is chosen on the left 
of 0, we take the point 1 above for a positive value of xi and 
be|ow for a negative one and the same for the other points. 
Mark a point above in the same distance from as P. Join 
the point P with the points 0, 1^ 2^ 3, 4, • • • and draw a broken 
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line OAoAiAiAzAa in such a manner that Aq is on the vertical 
through Ci and OAo is parallel to PO, Ax on the vertical through 
Cs and A^i parallel to PI, A% on the vertical through Cj and AiA% 
parallel to P2 and so on. Then the ordinate po of i^o will have 
the same length as Of^ and will be directed upward when the 
direction of oq is to the right, and downward when the direction 
of ao is to the left. The difference ^ -- po of the ordinates of Ai 
and Aq is equal in length to Oiari, as pi ~ Po &nd ai have the same 
ratio as 01 and PO. Ai will be above or below Aq according to 
the line aiXi being directed to the right or to the left and it is 
understood that aiZi has the same direction as ai for positive. 



//" 




dm OU^ 



Fig. 15. 

values of xi and a direction opposite to Oi for negative values 

of xi. Thus the ordinate yi has the same length as the line 

ao + aixi and its direction is upward or downward according to 

the direction of the line Oq + aiXi being to the right or to the left. 

In the same way it is shown that the ordinate y% of the point A% 

is equal in length to 

Oo + aiXi + OiXi, 
and 2^ to 

and so on, the direction upward or downward corresponding 
to the positive or negative value of the linear function. 
If the values o{ xi,QDi, • • • Xn satisfy the equation 

ao + aiXi+ 020^ + • • • + OnXn ■= 

the ordinate pn must vanish, that is to say, the point An must 
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coincide with C^Hf ^^ ^^^ o' the line an* And vice versa if ^» 
and Cn+i coincide the equation is satisfied. Consequently if we 
know all the values but one of the numbers Xi, qd^, '" Xn the 
unknown value can be found graphically. For suppose a:^ to be 
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the unknown value we can, beginning from 0, find the broken 
line as far as ^2 and beginning from the other end An we 
can find it as far as ils (Fig. 15). A parallel to A%A$ through P 
furnishes the point 3 on the axis of ordinates. li xuXt, • • • x».i 
are known and only Xn not, we can draw the broken line as far 
as An~i and as An has to coincide with Cn+i> we can draw a parallel 
to An^iAn through P and find the point n on the axis of ordinates 
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that determines the value Xn by the ratio On/PO or On/Oo. In 
Figs. 15 and 16 all the coefficients oo, ctu * * '9 are positive. A 
negative coefficient as is shown in Fig. 17. The only difference 
is that Ce lies to the left of C5 and consequently the broken line 
passes from Aa back to A^. 

If we keep the points 0, 1, 2, • • •, in their positions but change 
the position of P to P' (Fig. 18) and repeat the construction of 
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the broken line, we obtain OAf!A\A% • • • instead of OAf^\A% 
The ordinate yj^ of the point AJ[ is evidently 



00 
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y« =* ^ Dvi + ^1 DVi + • • • + fl« 
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and therefore 



ro 



PO 



FO 



P« ~ pfQ ?«• 



That is to say, by changing the position of P without changing 
the position of the points 0, 1, 2, • • • we can change the scale of 
the ordinates of the broken line. They change inversely pro- 




«oCi «i C;^(7,«i(74 ^*Q 



Fig. 18. 

portional to PO. It may be convenient to make use of this 
device in order to make the ordinates a convenient size inde- 
pendent of the scale that we have chosen for the points 0, 1, 2, • • • 
that determine the values 

01 02 

^i^OQ' ^=00'"'- 

A linear equation with only one unknown quantity 

oo + aiXi = 

is solved by drawing a parallel to AoAi through P. Let a second 
equation be given with two unknown quantities 

bo + biXi + (m = 0. 

The lines bo, bi, b^ Bre laid off as before. Knowing xi as the 
solution of the first equation we can construct the broken line 
OBoBi corresponding to the second equation and as B2 must 
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coincide with the end of 12, we can draw a parallel to B1B2 
through P and find a^. In a similar manner we can find xs 
from a third equation 

and so we can find any number of unknown quantities, if 
each equation contains one unknown quantity more than those 
before. 

In the general case when n unknown quantities are to be 
determined frogi n linear equations each equation will contain 
all the unknown quantities, and therefore we cannot find them 
one after the other as in the case just treated. But it can be 
shown that by means of very simple constructions the general case 
is reduced to a set of equations, such as has just been treated. 

^ ^^ ^^ Let us begin with two 

T oi V vtti ; 1 1 equations and two un- 

known quantities. 

ho + biXi + 622^ = 0. 




60 



V 

Fig. 19. 



The lines Oq, ai, 02 are laid off on a horizontal line OA0A1A2 and 
the lines feo> &i> h on another horizontal line 0'BoBiB2 (Fig. 19). 
Now let us join and 0', Ao and Bo, Ai and Bi, A2 and B2 by 
straight lines and let us draw a third horizontal line intersecting 
them in the points O'^CoCiCi. These points correspond to a 
certain linear function 

Co + CiXi + C2X2, 

and it can be shown that it vanishes when xi and QD2 are the same 
values for which the first two linear functions vanish. Let the 
distance of the first two horizontal lines be I and the distance of 
the third from the first and second h and k. Then it can readily 
be seen that 

Co = Oo "t" "T (&0 — fl^o) = T Clo + "T &0- 
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For a parallel to 00' through Aq defines with the line AoBq on 
the third and second horizontal line segments equal to Co -- ao 
and &o ~~ oo and as these segments have the ratio h/l, it follows 
that 

Co = oo + 7 (60 "- oo) = 7 oo + 7 60. 

By drawing a parallel to AoBq through Ai and to AiBi through 
At or through B2 (which comes to the same thing), we convince 
ourselves in the same way that 

Ci = ai + 7 (61 — ai) = 7 ^^ "t- 7 i| 
and 

Multiplying the equation 

Oo + (h^ + (h^ ™ 
by k/l and the equation 

bo + biXi + b^-O 
by h/l and adding the two products, we obtain 

Co + CiXi + Oi?D% = 0. 

The third horizontal need not lie between the first two. If it 
lies below the second we have merely to give k a negative value 
and if it lies above the first we have to give h a negative value 
and the same formulae for cq, Ci, cz hold good. Consequently the 
conclusion remains valid, that from the first two equations the 
third follows. 

Now as we are perfectly at liberty to draw the third horizontal 
line where we please, we can let it run through the intersection 
of the straight lines AiBi and A^Bt. In this case the points 0% 
and Ct must coincide and consequently Oz must vanish. If Ci 
does not vanish we can by what has been shown above find Xi 
and with Xi we can find qdz from either of the two first horizontal 
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lines. In case ci also vanishes, that is to say, in case the three 

straight lines A^B^, AiBi, AoBo all pass through the same point, 

while 00' does not pass through it, the two given equations 

cannot simultaneously be satisfied. For if they were, it would 

follow that 

Co + ciaji + (Wi = 0, 

and as ci and Oz are zero cq would have to be zero, which it is not 
as 00' is supposed not to pass through the intersection of AtBt, 
A\Bi and Af^Ba. If on the other hand all four lines AtB^^ AiBu 
AoBof 00' pass through the same point, cq, ci and c^ will all three 
vanish. In this case the two given equations do not contradict 
one another, but bt/hh will be proportional to 000102. The 

4r 
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second equation will therefore contain the same relation between 
Xi and QDi as the first, so that there is only one condition for Xi 
and 0^ to be satisfied. We- may then assign any arbitrary value 
to one of them and determine the value of the other to satisfy the 
equation. 

In the case of two linear equations of any number of quantities 
Xi, Xi, ••' XnVfe can by the same graphical method eliminate one 
of the quantities. In Fig. 20 this is shown for two linear equa- 
tions with six unknown quantities. The tWQ horizontal lines 
OAoAiAzAiA^AsAB and O'BJBiBzBzB^JBhBt represent two linear 
equations. Through. the intersection of AzBt and AJS^ a third 
horizontal line is drawn intersecting the lines 00', -4oBo, AiB^ 
• • • AfPt in 0"CqCi • • • Cg. As Cs and Ca coincide, the line C4 
vanishes and x^ is eliminated, so that the equation assumes the 

form 

Co + CiXi + CzX2+c^z + CtPCh + Ce^ft = 0. 
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Suppose now that a set of six equations with six unknown quan- 
tities is represented geometrically on six horizontal lines. We shall 
keep one of these; but instead of the other five we coiistruct five 
new ones from which one of the unknown quantities has been 
eliminated by means of the first equation. Now it may happen 
that at the same time another unknown quantity is eliminated, 
then this quantity remains arbitrary. Of the five new equations 
we again keep one that contains another unknown quantity and 
replace the four others again by four new ones from which this 
unknown quantity has been eliminated. Going on in this 
manner the general rule will be that with each step only one 
quantity is eliminated, so that at last one equation with one un- 
known quantity remains. Instead of the given six equations 
with six unknown quantities each, we now have one with six, 
one with five and so on down to one with one. The geometrical 
construction shows that this system is equivalent to the given 
system, for we can just as well pass back again to the given 
system. We have seen above how the unknown quantities 
may now be found geometrically. It may however happen in 
si)ecial cases that with the elimination of one unknown quantity 
another is eliminated at the same time. To this we may then 
assign an arbitrary value without interfering with the possibility 
of the solution. Finally all unknown quantities may be elimi- 
nated from an equation. If in this case there remains a term 
different from zero it shows that it is impossible to satisfy the 
given equations simultaneously. If no term remains, the two 
equations from which the elimination takes its origin contain the 
same relation between the unknown quantities and one of them 
may be ignored. 

§ 5. The Orapkical Handling of Complex Numbers. — A complex 
number 

z = «+ yi 

is represented graphically by a point Z whose rectangular coordi- 
nates correspond to the numbers z and y. The units by which 
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the coordinates are measured^ we assume to be of equal length. 
We might also say that a complex number is nothing but an 
algebraical form of writing down the coordinates of a point in a 
plane. And the calculations with complex numbers stand for 
certain geometrical operations mth the points which correspond 
to them. 
By the ''sum'' of two complex numbers 

3i = a?i + yii and Zi = x^ + yd 
we understand the complex number 

fs$^ xg + yd 
where 

xg- xi + xst and »§ = »i + ai, 
and we write 

2!| = Zl + %. 

Graphically we obtain the point Zz representing z^ from the 
points Zi and Z^ representing Zi and 2^ by drawing a parallel 

to OZ2 through Zi and making ZiP 
(Fig. 21) equal to 0Z% in length 
and direction or by drawing a paral- 
lel through Z2 and making 2^P 
equal to OZi in length and direc- 
tion. The coordinates of P are 
evidently equal to Xi + Xx and 

yi + Vi^ 
Two complex numbers z and z' 

Fxo. 21. are called opposite, when their sum 

is zero. 

z + z' = or x^ -- x' and y^^}f or z ^ -- 2^. 

The corresponding points Z and Z' are at the same distance from 
the origin but in opposite directions. 

The difference of two complex numbers is that complex 
number, which added to the subtrahend gives the minuend 

2f2 + (Zl — 2??) = 2l- 
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Therefore 

2i — % = (ici — a^) + (yi — v»)i. 

This may also be written 

zi + Z2' where %'= — 2^ = — o^ — jfei. 

That is to say^the subtraction of the complex number 2^ from zi 
may be effected by adding the opposite number — 2^. For the 
geometrical construction of the point Z corresponding to zi — 2^ 
we have to draw a parallel to OZ2 through Zi and from Zi in 
the direction from Z2 to we have to lay oflP the distance Z2O. 
Or we may also draw from a line equal in direction and in length 
to ZtZi. This will also lead to the point Z representing the 
difference zi — Zg. 

The rules for multiplication and division of complex numbers 
are best stated by introducing polar coordinates. Let r be the 
positive number measuring the distance OZ in the same unit 
of length in which z and y measure the abscissa and ordinate, so 

that 

r^ s/a? + f 

and let ^ be the angle between OZ and the axis of x, counted in 
the direction from the positive axis of x toward the positive 
axis of y through the entire 
circumference (Fig. 22). Then 
we have 

a;=rcos^, y=rsin^ 

and 

z— ar+yi = f(cos^+sinv>0.. 

Let us call r the modulus 
and ^ the angle of z. The an- yiq. 22. 

gle may be increased or di- 
minished by any multiple of four right angles without altering 
z, but any alteration of r necessarily implies an alteration of z. 
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Aoooiding to McMVie's theorem, we can write 

By the product of two complez numbers 

zi = fi^* and a^ = v**^ 

we understand that complez number zg whose modulus ri is 
equal to the product oi the moduli fi and fs and whose angle ^ 
is the sum of the anj^es ^ and ^ or differs from the sum only by 
a multiple of four right anj^es 

z^ — ZA = fiiiic^**"*"*^. 

The definition of division follows from that of multiplica- 
tion. The quotient zi divided by :^ is that complex number, 
which multiplied by 2^ gives zi. Therefore the product of its 
modulus with the modulus of 2^ must be equal to the modulus of 
Zi and the sum of its angle with the angle of z% must be equal to 
the angle of zi. Or we may also say the modulus of the quotient 
zi/zi is equal to the quotient of the moduli fi/rs and its angle is 
equal to the difference of the angles ^ — ^. An addition or 
subtraction of a multiple of four right angles we shall leave out 
of consideration as it does not affect the complex number nor 
the point representing it. 

The geometrical construction corresponding to the multi- 
plication and division of complex numbers is best described by 
considering two quotients each of two complex numbers that 
give the same result. Let us write 

zi/zi = aii/zi. 

The geometrical meaning of this is that 

ri/fa = r8/r4, 
and 

<pi — <ft = ^ — ^ 
That is to say, the triangks Z1OZ2 and Z^OZi aie geometrically 
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similar (Fig. 23). When three of the points Zi^ Z^, Zz, Zi are 
given the fourth can evidently be found. For instance let 
Zu Zt, Z4 be given. Draw a parallel to Z1Z2 intersecting 02s 
at a distance u from 0. This point together with the inter- 
section on OZi and with will form the three comers of a tri- 
angle congruent to the triangle Z^^iO. It will be brought into 





Fig. 23. 
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the position of Zi^zO by being tiumed round so as to bring the 
direction of the side in 0Z% into the position of OZ4. Thus the 
direction of OZt and its length may be found. 

This construction contains multiplication as well as division as 
special cases. Let Z4 coincide with the point x = 1^ p = 0, so 
that 24 = 1 (Fig. 24), then we have 

Zilz^ ^ 2(j or 2i = 2^. 

From any two of the points Zi, Z2, Zz a simple construction ^ves 
us the third. 

The geometrical representation of complex numbers may be used 
to advantage to show the properties of harmonic oscillations. 

Let a point P move on the axis of x^ so that its abscissa at the 
time t is given by the formula 

a: = r cos (n/ + a), 

n, r and a being constants. We call r the amplitude and nf + a 
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the phase of the motion. The point P moves backwards and 
forwards between the limits x^ r and x^ ^ r. The time 
T ^ 2r/n is called the period of the oscillation, it is the time in 
which one complete osdllation backwaids and forwards b per- 
formed. 
Now instead of x let us consider the complex number 

a = r cos(n/ + a) + r sin(n/ + a)i 
or 

of which X is the abscissa and let us follow the movement of 
the point Z. For t^ Owe have 

Deagnating this value by Zo, we can write 

2 = 20^***. 

The geometrical meaning of the product 

is that the line OZq is turned round through the angle tU. For 
the modulus of e*** being equal to 1 the modulus of zq is not 

changed by the multiplication. The 
movement of the point Z therefore 
consists in a uniform revolution of 
OZ round 0. At the moment <= 
the position is OZq and after the 
time T = 2x/n the same position is 
occupied again. The revolution goes 
on in the direction from the positive 
axis of a; to the positive axis of y 
(Fig. 25). 

The movement of Z is evidentiy 
simpler than the movement of the 
projection P of Z on the axis of x. 
Let us consider a motion composed of the sum of two harmonic 







Fig. 26. 
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motions of the same period but of different amplitudes and 

phases 

X ^ ficos (fd + cti) + fgC08(fd + as), 

and let us again substitute the motion of the point Z correspond- 
ing to the complex number 

For < » the first term is 
and the second term 

Introducing Zi and sk into the expression for z we have 

z = aie"** + %c*^ «= (3i + zt)e*^ = V*** 
where 

2|B = 2fl + a^. 

This sho¥rs at once that the movement of Z is a uniform circular 

movement consisting in a uniform revolution of OZ round 0. 

The position at the moment < = is OZz corresponding to the 

complex number 

2iB = «i + asi. 

The projection of Z on the axis of z has the abscissa 

where ff and 0$ designate modulus and angle of 2iB. Thus the 
sum of two harmonic motions of the same period is shown also 
to form a harmonic motion. 

The same holds for a sum of any number of harmonic motions 
of the same period. For the complex number 

where ri, fi, •*• fx; cti, a^, -•' ax and n are constants may be 

written 

z » zie^ + 8^*^ + • • • + M"** 
or 
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where 

«• = »i + ak + h 2x. 

The movement of Z therefore, exeeptmg the case 80 »= 0, oonasts 
in a uniform revolution of OZ round 0, OZ always keeping the 
same length equal to the modulus of Zo* The position of OZ at 
the moment < = is OZ^ 
The motion of a point P whose abscissa is 

X s= ojT^ cos (n< + a) 

where a, Ic^ n, a are constants (a and k positive) is called a damped 
harmonic motion. It may be looked upon as a harmonic motion, 
whose amplitude is decreasing. To study this motion let us 
again substitute a complex number 

z = (UT^ cos (fd+ a) + a^ sin (n< + ct)i, 
or 

or 



where zo is written for the complex constant o^. 
The product 

is a complex number corresponding to a point Zi on the same 
radius as Zq, coincident with Zo at the moment i ^ but ap- 
proaching in a geometrical ratio after t=^0. In unit of time 
the distance of Zi from decreases in the constant ratio ^ :1. 
The multiplication with e*^ turns OZi round through an angle 
ni. We may therefore describe the motion of Z as a imiform 
revolution of OZ round 0, Z at the same time approaching 
at a rate uniform in this sense that in equal times the distance 
is reduced in equal proportions (Fig. 26). At the moment < = 
the position coincides with Zq. We speak of a period of this 
motion meaning the time T = 2T/n in which OZ performs an 
entire revolution round 0, although it does not come back to its 
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original position. Any part of the spiral curve described by Z 
in a given time b geometrically similar to any other part of the 
curve described in an interval of equal duration. For suppose 
the second interval of. time hap- 
pens r units of time later, we 
shall have for the first interval 




and for the second interval 

Now if zi and 2^ are the values 
of z at two moments ti and U of 
the first interval and zi and zg' 
the corresponding values of z' 
at the moments ii+ r and ^ + r of the second interval, we have 



1 
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^ ss ^(«i-li).^»(fa-«i)< = ?L 
Zi Zi'' 

Therefore the triangle Z1OZ2 is geometrically similar to the 
triangle ZiOZ^'. As Zi and Z% may coincide with any points 
of the first part of the curve, the two parts are evidently geo- 
metrically similar. 

The projection of Z on the axis of x performs oscillations 
decreasing in amplitude. The turning-points correspond to those 
pK>ints of the spiral curve described by Z, where its tangent is 
parallel to the axis of y, that is to say, where the abscissa of dzjdt 
vanishes. 



Now 



dz 

•t: = »o(— i + m)€^e'^ = (— i + m)z 



or 



dz 

z 



= — i + tii = ptf^ 
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where p and X are the modulus and angle of the complex number 

Consequently, if we represent dzjdi by a point Z', the triangle 

Z'OZ will remain geometrically similar to itself. The turning 

points of the damped oscillations correspond to the moments 

when OZ' is directed vertically upward or downward or when the 

angle of dzjdi is equal to x/2 or 3t/2. The angle of z will then 

be t/2 — X or 3x/2 — X plus or minus any multiple of 2t. As 

the angle of z, on the other hand, b changing in time according 

to the formula 

fd + a, 

we find the moments where the movement turns by the equation 

ni + a^ tI2-\+2Nt, 
or 

fd + a^ 3t/2 - X + 2Nt, 

N denoting any positive or negative integral number. The time 
between two consecutive turnings is therefore equal to w/n, that 
isy equal to half a period. All the points Z corresponding to 
turning points lie on the same straight line through the origin O 
forming an angle 3x/2 — X with the direction of the positive axis 
of X. The amplitudes of the consecutive oscillations therefore 
decrease in the same proportion as the modulus of z, that is 

to say, in half a period in the ratio ^ir* 

Let us consider the vibrations of a system possessing one 
degree of freedom when the system is subjected to a force varying 
as a harmonic function of the time and let us limit our considera- 
tions to positions in the immediate neighborhood of a position 
of stable equihbrium. If the quantity x determines the position 
of the system the osdllations satisfy a differential equation of 
the form 

where m, i, n, p, F are positive constants. 

^ See for instance Rayleigh, Theory of Sound, Vol. I, chap, in, { 46. 
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This is another case where the introduction of a complex 

variable 

z = x + yi 

and the geometrical representation of complex numbers helps to 
form the solution and to survey the variety of phenomena that 
may be produced. 

In order to introduce 2 kt us simultaneously consider the 
differential equation 

m^ + iJ+nV =Fsinp<, 

and let us multiply the second equation by i and add it to the 
first. We then have 

The movement of the point Z representing the complex number 
z then serves as well to show the movement corresponding to x. 
We need only consider the projection of Z on the axis of x. 

A solution of the differential equation may be obtained by 
writing 

Introducing this expression for z and cancelling the factor e^** 

we have 

Zo(— mp^ + kpi + n*) == F, 
or 

^ F , 

zo is a complex constant, that may be represented geometrically 
as we shall see later on. 
This solution 

Z = ZoC^" 

is not general. If z' denotes any other solution so that 
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we find by subtracting the two equations 

m '^ ' + fe ' ^ ' + nV - a) = 

or writing 

a' — 2 =s ti, 

d*M _ du ^ ^ 

^W^ ~di'^ 
The general solution of this equation is 

where t^i and i/s are arbitrary constants and Xi and Xs are the 
roots of the equation for X 

mX* + ifcX + n* = 0, 



X«J" 2m \4m« ^" 



If y/4m ' is greater than n^ so that the square root has a real 
value^ Vk^l^rn^ — n* will certainly be smaller than k/2m. There- 
fore Xi and Xe will both be negative and the moduli of the complex 
numbers ui^^ and ihfi^ will in time become insignificant. If, 
on the other hand, l^fArn^ is smaller than n^ both complex 
numbers Uie^* and u^^ correspond to points describing spirals that 
approach the origin, as we have seen above, in a constant ratio 
for equal intervals of time. Therefore they will also in time 
become insignificant. 
After a certain lapse of time the expression 

will therefore suffice to represent the solution. 

The point Z moves uniformly in a circle round of a radius 
equal to the modulus of zq, completing one revolution in the 
period 2x/p, the period of the force acting on the system. The 
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movement of the projection of Z on the axis of x is given by 

ar = fo cos (jrf + a), 

where fo is the modulus and a the angle of zo. It is a harmonic 
movement with the same period as that of the force* Fcosp^, 
but with a certain difference of phase and a certain amplitude 
depending on the values of F, m, i, n, p. 

It is important to study this relation in order to survey the 
phenomena that may be produced. For this purpose the geo- 
metrical representation of complex numbers readily lends itself. 

In the expression for Zq 

let us consider the denominator 

and let us suppose the period of the force acting on the system 
not determined, while the constants of the system m, h^ n and 
the amplitude of the force F have given values. The quantity p 
is the number of oscillations of the force during an interval of 
2x units of time. This quantity p we suppose to be indeter- 
minate and we intend to show how the amplitude and phase 
of the forced vibrations compare with the amplitude and phase 
of the force for different values of p. 

Let us plot the curve of the points corresponding to the complex 
number 

n* — mi^ + ipi, 

where p assumes the values p = to +00. 

This curve is a parabola whose axis coincides with the axis of 
X and whose vertex is in the point a: = n^ y = 0. We find its 
equation by eliminating p from the equations 

a; = n* — mp*, y = ftp, 
viz.. 
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« = n» — j^jr. 

But it is better not to eliminate p and to plot the different points 
for different values of p. In Kg. 27 the curve is drawn f or p = 

to 3 and the points for 
p==0, 1, 2, 3 are marked. 
The ordinates increase in 
proportion to p; they are 
equal to 0, k, 2k, 3k for 
p = 0, 1, 2, 3. The dis- 
tance between the projec- 
tion of any point of the 
curve on the axis of x and 
It is equal to 0, m, 4m, 9m for 




Feo. 27. 



the vertex is proportional to p*. 
p = 0, 1, 2, 3. 

For any point P on the parabola let us denote the distance 
from by r and the angle between OP and the positive axis of 
;r by ^ so that 



Then we have 
and consequently 

and 



F 
r 



« = ~ cos (p^ — ^). 



The amplitude F/r of the forced vibration is inversely propor- 
tional to r. Thus our Fig. 27 shows us what the period of the 
force must be to make the forced vibrations as large as possible. 
It corresponds to the point on the parabola whose distance from 
is smallest. It is the point where a circle round touches the 
parabola. In Fig. 27 this point is marked JR. It may be called 
the point of majdmum resonance. When the constants of the 
system are such that the ordinate of the point, where the parabola 
intersects the axis of y is small in comparison with the abscissa 




GRAPHICAL CALCULATION. 39 

of the vertex, then OR will lie close to the axis of y (Fig. 28). In 
this case the angle between OR and the positive ans of x will be 
very nearly equal to 90^, that is to say, the forced oscillations will 
lag behind the force oscil- 
lations by a little less than a f ^ 
quarter of a period. Keep- 
ing m and n constant, this 
will take place for small val- 
ues of kf t. e., for a small 
damping influence. A small j^^^ 28. 
deviation of p from the fre- 
quency of maximum resonance will throw the point P away from R, 
so that r increases considerably and ^ becomes either very small 
(for values of p smaller than the frequency of maximum resonance) 
or nearly equal to 180** (for values of p larger than the frequency 
of maximum resonance). In other words for small values of k the 
maximum of resonance is very sharp. A deviation of the period 
of the force from the period of maximum resonance will lessen the 
amplitude of the forced vibration considerably. The lag of its 
phase behind that of the force will at the same time nearly vanish, 
when the frequency of the force is decreased or it will become nearly 
as large as half a period, when the frequency of the force is in- 
creased. For larger valuea of k the parabola opens out and this 
phenomenon becomes less marked. The minimum of the radius r 
becomes less pronounced. The angle between OR and the axis of 
X becomes smaller and smaller and for a certain value of k and all 
larger values the point R will coincide with the vertex of the para- 
bola. In this case, there is no resonance. When the period of 
the force increases indefinitely (p becoming smaller and smaller) 
the amplitude of the forced vibration will increase and will 
approach more and more to the limit 

F 

but there will be no definite period for which the forced vibra- 
tions are stronger than for all others. 
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The Graphical Representation of Functions of One ob 

More Independent Variables. 

§ 6. Functions qf One Independent Variable. — ^A function y of 
one variable x 

» = /(«; 

b usuaUy represented geometrically by a curve, in such a way 
that the rectangular coordinates of its points measured in certain 
chosen units of length are equal to x and y. This graphical rep- 
resentation of a function is exceedingly valuable. But there is 
another way not less valuable for certain purposes, more used in 
applied than in theoretical mathematics, which here will occupy 
our attention. 

Suppose the values of y are calculated for certain equidistant 
values of x, for instance: 

ar = - 6, - 5, - 4, - 3, -2, - 1, 0, 

+ l, + 2, + 3, + 4, + 5, + 6, 

and let us plot these values of y in a uniform scale on a straight 

line. Draw the uniform scale on 
one side of the straight line and 
mark the points that correspond 
to the calculated values of y on 
the other side of the straight line. 
Denote them by the numbers x 
that belong to them (Fig. 29). 
The drawing will then allow us to 
read off the value of y for any of 
the values of x with a certain ac- 
curacy depending on the size of the 

scale and the number of its partitions and naturally on the fine- 

40 



7- 

•--4 

"y^i — 



Feo. 20. 



GHAPmCAL REPKE8ENTATI0N OF FUNCTIONS. 



41 



ness of the drawing. It will also allow us to read off the value 
of y for a value of x between those that have been marked, if 
the intervals between two consecutive values of x are so small 
that the corresponding intervals of y are nearly equal. We can 
with a certain accuracy interpolate values of x by sight. On the 
other hand, we can also read off the values of x for any of the 
values of y. We shall call this the representation of a function 
by a scale. 

We can easily pass over to the representation of the same 
function by a curve. We need only draw lines perpendicular 
to the line carrying the scales through the points marked with 
the values of x and make their length measured in any given 
unit equal to the numbers x that correspond to them (Fig. 29). 

In a similar way we can pass 
from the representation of the 
function by a curve to the rep- 
resentation by a scale. 

The representation by a scale 
may be imagined to signify the 
movement of a point on a straight 
line, the values of x meaning the 
time and the points marked with 
these values being the positions 
of the moving point at the times 
marked. By passing over to the curve tbe movement in the 
straight line b drawn out into a curve with the time as abscissa 
(Fig. 30). 

The representation by a scale is used in connection with the 
representation by a curve for the purpose of drawing a function 
of a function. 

Let y be a function of x and x a function of t. Then we wish 
to represent y as a function of L 

Let y = Six) be given by a curve in the usual way and let 
X = ip{i) be given by a scale on the axis of x marking the points 
where < = 0, 1, 2, • • •, 12. We then find the values of y corre- 
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sponding to the values < = 0, 1, 2, • • • , 12 by drawing the oidi- 
nates of the curve y = S(x) for the abscissas marked < = 0, 1, 2, 




V=r/(**(t» 



Fio. 31. 



• • -1 12. These ordinates as a rule will not be equidistant But 
as soon as we move them so as to make them equidistant, they 

form the ordinates of the curve 

V = /(v(0) 
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1 



with i as abscissa (Fig. 31). 

The representation of a func- 
tion by a scale may be general- 
ized in the respect that ndther 
of the two scales facing one an- 
other on the straight line need 
necessarily be uniform. The in- 
tervals of both scales may vary 
from one side of the scale to the 
other. If the variation is suflS- 
ciently slow the interpolation can nevertheless be effected with 

accuracy. We may look at this case as composed of two cases 
of the first kind. 

fix) = y and y = f(0- 
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These scales are placed together, so that the scale x touches the 
scaled 

while the scale y is cut out (Fig. 32). 

§ 7. The Principle of the Slide RtUe. — ^Let us investigate how 
the relation between x and t changes by sliding the x- and ^scales 
along one another. 

If we slide the ovscale through an amount y = c so that a 

point of the ovscale that was opposite to a certain point y of the 

^-scale, now is opposite y + c, then the relation between x and t 

represented by the new position of the scales will be given by 

the equation 

fix) = g(t) + c. 

II X, t and s^, ff denote two pairs of values that are placed 
opposite to one another, we shall have simultaneously 

fix) = ^(0 + c, 

fix')^gif) + c, 
or by eliminating c 

fix) - git) = fix') - gif). 

The ordinary slide rule carries two identical scales y = Ipg a; and 
2^ = log ^ that are able to slide along one another, x and t running 
through the values 1 to 100. We therefore have 

log a: — log < = log x' — log f, 
or 

t^ f' 



I A i* f V 
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That is to say, in any position of the x^ and ^scale any two values 
X and i opposite each other have the same ratio (Fig. 33). This 
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is the princiide oa which the use 6! the sfide rule b founded. 
It enables us to calcuUte any of the four quantities x, i, x', t' 
if the other three are given. Suppose, for example, x, t, x' 
known. We set the scales so that x appears opposite to i. 
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then ^ is read off opposite to x'. On the other edges the slide 
rule carries two ^milar scales one double the size of the other 
(Fig. 34). We may write 

y = 2 log Z and y « 2 log T. 

By means of a little frame carrying a crossline and sliding over 
the instrument, we can bring the scales x and T ori and X op- 
posite each other. If, for example, for any position of the 
instrument x, T and x', T are two pairs of values opposite each 

other, then 

logx - 21og T = logar' - 21og T\ 
or 

If any three of the four quantities x, T, x', T are known the 
fourth may be read off. Thus we find the value 

xT^ 



by setting T opposite to x and reading off the value opposite to 
T\ Or we can find the value of 



4 



X 



by setting x opposite to T and reading off the value opposite oc^. 
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Let US reverse the part that carries the scales t, T so that x 
slides along T and X along t, but in the opposite order (Fig. 35). 




Fig. 35. 

The scales t, T may then be expressed by 

« 

y^l-logt and y = / - 2 log T, 

/ being the entire length of the scales. 

By setting the instrument to any position and considering the 
scales X and / or X and T by means of the cross line we have 
loga;+ log* = loga/ + loe<' and logZ+ log T = logZ' + log T' 

or 

xt-^x't' and XT^^rr, 

so that any two values opposite to one another have the same 
product. 

For X and T we have 

loga; + 21og 7= log a:' + 2 log T, 
or 

xT = x'r\ 

Let us apply this to find the root of an equation of the form 

u^ + au ^ b. 
Divide by t^ so that 

tt' + a = - 
u 

and set r = 1 opposite to X » 6. Then taking 7^ = t^ we find 
on the same cross line t^u^ and X = b/u, so that we read the 
two values u^ and b/u directly opposite to each other on the 
scales t and X. If b/u is positive, it decreases while u^ increases. 
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Running our eye along we have to find the place where the differ- 
ence hju — u^ is equal to a. Having found it the T-^cale gives 
us the root of the equation. For example take 

w» - 5tt = 3, 
or 

u 

We set r = 1 opposite X = 3 and run our eye along the scales 
X and t (Fig. 36), to find the place where ^ — 5 = X. We find 
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Fig. 30. 

it approximately at ^ = 6.2, and on the 7^--scale we read off 
T = 2.50 as the approximate value of the root. This is the 
only positive root. But for a negative root 3/u is negative, 
and therefore the positive value of 3/m plus u^ would have to be 
equal to 5. We run our eye along and find t = 3.37 opposite to 
X = 1.63, approximately corresponding to T = 1.84. There- 
fore — 1.84 is another root. As the coefficient of u^ in the first 
form of the equations vanishes it follows that the sum of the 
three roots must be equal to zero. This demands a second 
negative root approximately equal to — 0.66. To make sure 
that it is so, we set the instnunent back and take the other end 
of the T-scale as representing the value T = 1 and give it the 
position this end had before. Running our eye along the 
scales X and t, we find t = 0.43 opposite to X = 4.57, giving 
Z + ^ = 5.00. On the T-scale we find 0.655, so that the third 
root is found equal to — 0.655. 

When b is negative there is always one and only one negative 
root. For u running through the values w = to — oo, u^-^b/u 
will run from — oo to + oo without turning. When b is positive 
there is always one and only one positive root; for then u^ — b/u 
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runs from — oo to + oo for w = to + oo. In the first case 

there may be two positive roots or none; in the second case there 

may be two negative roots or none. For positive values of a 

one root only exists in either case. This is easily seen in the first 

form of the equation 

if + au = b, 

because from a positive value of a it follows that u^ + au will 
for M = — 00 to + 00, run from — oo to + oo without turning 
and will therefore pass any given value once only. 

In order to decide whether in the case of a negative value of a 
there are three roots or only one let us write 

w* — - = — a. 
u 

For negative values of b we have to investigate whether there 
are positive roots. For positive values of u the function u^^b/u 
has a minimum, when the differential coefficient vanishes, i. e., for 

2w + -i = 0, 
or 

u 

Having set our slide rule so that t gives us u^ and X gives us 
— b/Uy we find the value u where the minimum takes place by 
running our eye along and looking for the values X, t opposite 
each other for which X is twice the value ott 

Then t+ X is the minimum of u^ — b/u, so that there will be 
two or no positive roots according to t + X being smaller or 
larger than — a. For positive values of b, we have to find out 
whether there are negative roots. The criterion is the same. 
After having set T = 1 opposite to , b and having found the 
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positive root, we find the place where 

Then ^ + Z is the tninimiiTn of all values that i^ — b/u assumes 
for negative^ values of u. If the minimum is smaller than — a 
there are two negative roots; if it is larger there are none. If it 
is equal to — a the two negative roots coincide. 
For the equation 

u 

for instance, we find t = 1.31 opposite to X = 2.62 (Fig. 36), 
so that 2t = 2.62 = X. Now ^ + Z = 3.93 is smaller than 5, 
therefore li* — d/u will assume the value 5 for two negative 
values of u on either side of the value li = — T = — 1.143 
for which the minimum of u^ — S/u takes place. 

On the same principle as the slide rule many other instruments 
may be constructed for various calculations. In all these cases 
we have for any position of the instrument 

fix) - g{t) = fix') - gin. 

where x, t are any readings of the two scales opposite each other 
and x'f the readings at any other place. f(x) and g(t) may be 
any functions of x and t It will only be desirable that they 
be limited to intervals of x and t, which contain no turning 
points. Else the same point of the scale corresponds to more 
than one value of a; or ^ and that will prevent a rapid reading 
of the instrument. 

Let us design an instrument for the calculation of the increase 
of capital at compound interest at a percentage from 2 per cent, 
upward. If x is the number of per cent, and t the number of 
years, the increase of capital at compound interest is in the pro- 
portion 

+ m)'- 
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We can evidently build an instrument for which 



0+1^)'= (i+0- 



For taldng first the logarithm and then the logarithm of the 
logarithm, we obtain 

l0g< + l0gl0g(l + ^) = lQg<'+lQglQg(l+ ~j. 

We have only to make the x-scale 

y = +loglog(l+ j^)-logIog(l+j|). 

and the ^scale 

y = log n — log t 

For a; = 2 we have y = and therefore in the normal position 
of the instrument t ^ n. On the other end we have ^ = 1 and 
therefore y = log n. Now let us take n = 100, so that y = 2 
for t == 1. Say the length of the instrument is to be about 24 
cm., then the unit of length for the y-scale would have to be 12 
cm. In the normal position of the instrument the readings x, t 
opposite to each other satisfy the equation 



0+roo)'=0+iio)" 



Opposite t = l,we read the value zi = 624 and this gives us 



( 



2 \^^ r 

l+ioo) =1 + 1^=1 + 6.24 = 7.24. 



A capital will increase in 100 years at two per cent, compound 
interest in the proportion 7.24 : 1. Or we may also say the 
number xi = 624 read off opposite ^ = 1 is the amount which is 
added to a capital equal to 100 by double interest of 2 per cent, 
in 100 years. The same position of the instrument gives us the 
number of years that are wanted for the same increase of capital 
5 
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at a higher percentage. For all the values x, t opposite to each 
other satisfy the equation 



(i + ^§5) = 7^ 



For any other given percentage x i^id any other given number 
of years t the increase of capital is found by setting x opposite 
to i and reading the x-«cale opposite to ^ = 1. The only restric- 
tion is that the ratio is not greater than 7.24, else ^ = 1 will 
lie beyond the end of the a?-scale. 

For a given increase of capital the instrument will enable us 
either to find the number of years if the percentage is given, or 
the percentage if the number of years b given, subject only to 
the restriction mentioned. 

We can build our instrument so as to include greater increases 
of capital by choosing a larger value of n. n == 1000, for in- 
stance, will make y = 3 for / = 1. If the instrument is not to 
be increased in size the scales would have to be reduced in the 
proportion 2 : 3. 

Let us consider another instance 

1 1 1 

» = i' »=»-?• 

In the normal position of the instrument the scale division 
marked x = oo corresponds to y ^ and is opposite to ^ = n. 
If we have / = oo on the other end, the length of the instrument 
will correspond to y ^ 1/n. Let us choose n » 0.1, so that the 
length of the instrument is y = 10. That is to say, the imit of 
length of the ^scale is one tenth of the length of the instrument. 
For any position of the instnmient we have 

If the scale division marked a; = oo is opposite to ( =: c we can 
write x' == 00, t' = c and have 
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< g > < t i 



X i c 

The instrument will therefore enable us to read off any one of 
the three quantities z, t, c, if the other two are given, the only 
restriction being that all three lie within the limits 0.1 to oo. 
The instrument may be used to determine the combined resistance 
of two parallel electrical re- 
sistances, for the resistances 
satisfy the equation 

R Ri R^ 

^, . Fio. 37. 

Similarly it may be used 
to calculate the distances of an object and its image from the 
principal planes of any given system of lenses. For if / is the 
focal length and x and i the distances of the object and its im- 
age from the corresponding principal planes (Fig. 37), the equar 
fion is 

On the back side of the movable part of an ordinary slide rule 
there generally is a scale 

y == 2 + log sin t 

When this part is turned round and the scale is brought into 

contact with the scale 

y = log X, 

we obtain for any position of the instrument 

log a; — log sin / = log «' — log sin f, 



or 



a:' 



sin t sint" 



for any two pairs of values x, i that are opposite each other. 
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Given two sides of a triangle and the angle opposite the larger 
of the two the instrument gives at once the angle opposite the 
other side. Similarly when two angles and one side are given, 
it gives the length of the other side. 
If a:' =5 a is the value opposite to f = 90®, we have 

X = asm L 

Thus we can read the position of any harmonic motion for any 
value of the phase. 
An instrument carrying the scales 

y » log sin x and y = log sin t 

enables us to find any one of four angles x, t, a^, f for which 

sin x sin xf 



sin ^ sin ^' 

if the other three are given. Thus, knowing the declination, 
hour angle and height of a celestial body, we can read the azimuth 
on the instrument. We have only to take x = 90® — height, 
t = hour angle, x' = 90® — declination, then f = azioduth or 
180® - azimuth. 

It is not necessary to carry out the subtraction 90® — height and 
90® — declination. The difference may be counted on the scale 
by imagining 0® written in the place of 90®, 10® in the place of 
80® and so on and counting the partitions of the scale backwards 
instead of forward. 

§ 8. Rectangular CoSrdinates with Intervals qf Varying Size. — 
The two methods of representing the relation between two 
variables either by a curve connecting the coordinates or by 
scales facing each other lead to a combination of both. 

Suppose the rectangular coordinates x and y are functions of 

u and V, 

X = ^{u) and y = ^(t»). 

The function x = <p(u) is represented by a uniform scale for x 
on the axis of abscissee facing a non-uniform scale for u. The 
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function y = i^(v) is represented by a uniform scale for y on 
the axis of ordinates facing a non-uniform scale for v. Tlirough 
the scale-divisions u let us draw vertical lines, and through the 
scale-divisions n let us draw horizontal lines. These two systems 
of parallel lines form a network of rectangular meshes of various 
sizes (Pig. 3S), and any equation between u and v may be repre- 
sented by a curve in this plane. 

The usefulness of this method will be seen by some examples. 
It enables us by a clever choice of the functions ^(u) and V'(i') 




to Amplify the form of the curve. It is easily seen, for instance, 
that a curve representing an equation /(u, s) = may always be 
replaced by a straight line, if we choose the ti-sc&le properly. 
For when the points u = 1, 2, 3, 4, ■ - ■ of the curve are not on 
a stmght line, let them be moved to a straight line without 
altering their ordinates (Fig. 39). This will change the u-scale 
but it will not alter the equation /(«, v) = now represented by 
the straight line. 
Suppose we want to represent the relation 



where a and b are given numbers. If u and v were ordinary 
rectangular coordinates the curve would be an ellipse. But if we 
make 

I = u* and y = if 
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the equation ci the line in rectangular coordinates becomes 

and the carve will therefore be a straight line running from a 
point on the positive axb ct xto u point on the positive axis of 
y. The point on the axis of x corresponds to the value ii = =^ a 

on the upscale, and the 
point on the axb of y cor- 
responds to the value v = 
sfe b on the 9-scale (Fig. 40). 
Any point on the straight 
line corresponds to four 
combinations +u, +r;— tx, 
+ t; M, — «; — tt, — V, be- 
cause X has the same values 
for opposite values of u 
and y for opposite values of 9. We can read v as a function of 
t^ or li as a function of v. 
If a second equation 

— 4- — s= 1 

is given, we find the conmion solutions of the two equations by 
the intersection of the corresponding straight lines. Fig. 40 
shows the solutions of the two equations 




-t^odD 



Fia.4a 



and 



— + ~ = 1 
22 T" 32 A 

42 + 52 - 1» 



approximately equal to w = =*= 1.2 and « = =*= 2.4. 
Another function much used in mathematical physics 



1? = Oe mi 
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may also be represented by a straight line by means of the same 

device. 

By making 

y = log r, a: = w^ 
we obtain 

y = loga-^, 

where logt; and log a are the natural logarithms of v and a. 
The ti'-scale is laid off on the axis of x and the i>-scale on the axis 
of y and we have to join the points t^ = 0, t? = a and u = m, 

V = a/e. The point v = a/e is found by laying off the distance 

V = 1 to t» = « from V =^ a downward (Fig. 41). We are not 
obliged to take the same units of length for x and y. 










v^oe 



Fio. 41. 



Suppose we had to i^d the constants a and m from two equa- 
tions 

and 

ui = aer^' 

Our diagram would furnish two points corresponding to uu ^ 
and Vt, Vt. The straight line joining these two points intersects 
the axis of ordinates at t^ = a and intersects the parallel through 
V = afe to the axis of abscissae at u = m. 
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In applied mathematics the problem would as a rule present 
itself in such a form that more than two pairs of values u, v 
would be given but all of them affected with errors of observation. 
The way to proceed would then be to plot the corresponding 
points and to draw a straight line through the points as best we 
can. A black thread stretched over the drawing may be used to 
advantage to find a straight line passing as dose to the points 
as possible (Fig. 42). 

In several other cases the variables u and v are connected with 
the rectangular coordinates x and y by the functions 

X ^logu and y = log v. 




10 
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♦■ 
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Fio. 42. 



Fia. 43. 



''Logarithmic paper'' prepared with parallel lines for equidbtant 
values of u and lines perpendicular to these for equidistant values 
of t^ is manufactured commercially (Fig. 43). 
By this device diagrams representing the relation 

where r, 9, c are constants are given by straight lines. For by 
taking the logarithm we obtain 

rx + sy =^ log c. 

The straight line connects the point li = c^ on the upscale with 
the point « = c*^* on the »-scale. 
Logarithmic paper b further used to advantage in all those 
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cases where a variety of relations between the variables u and v 
are considered that differ only in u and v being changed in some 
constant proportion. If u and v were plotted as rectangular 
coordinates the curves representing the different relations be- 
tween u and V might all be generated from one of them by altering 
the scale of the abscissae and independently the scale of the ordi- 
nates, so that the appearance of all these curves would be very 

different. Let us write 

f{u, v) = 0, 

as the equation of one of the curves. The equations of all the 
rest may then be written 

/ (." B ) - »• 

where a, b are any positive constants. The points u, v of the 
first curve lead to the points on one of the other curves by taking 
u a times as great and t^ 6 times as great. For if we write ti' = aii 
and f/ == bu the equation f(u, v) =^ leads to the equation 
between u^ and v': 



/(^.f) = o. 



Using logarithmic paper the diagram of all these curves be- 
comes very much simpler. The equation /(w, t») = is equivalent 
to a certain equation <p(x, y) = 0, where a: = log w, y = log v. 
Now let x', y' be the rectangular coordinates corresponding to 
u't i/ so that 

xf = log ix' = log t^ + log a = a; + log a, 

yi = log t»' = log t» + log 6 = y + log 6. 

The point ar', y' is reached from the point ar, y by advancing 
through a fixed distance log a in the direction of the axis of x 
and a fixed distance log h in the direction of the axis of y. The 
whole curve 

M t^) = 
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drawn on logarithmic paper is therefore identical with aO the 
curves 



f{Vl)-o- 



It can be made to coincide with any one of the curves by 
moving it along the directions of x and y. 

§ 9, Functions cf Two Independent Variables. — ^When a func- 
tion of one variable y = f(x) is represented by a curve, the values 
of X are laid off on the axis of x and the values of y are represented 
by lines perpendicular to the axis of x. In a similar way a 
function of two independent variables 

2 = /(«> y) 

may be represented by plotting x and y as rectangular coordinates 
and erecting lines perpendicular to the xy plane, in all the 
points X, y, where f{x, y) is defined and making the lengths of 
the perpendiculars proportional to z. In this way the function 
corresponds to a surface in space. Now there are practical 
difficulties in working with surfaces in space and therefore it 
appears desirable to use other methods, that enable us to represent 
functions of two independent variables on a plane. This may 
be done in the following way. 

Taking x, y as rectangular coordinates all the points for which 
/(x, y) has the same value form a curve in the xy plane. Let 
us suppose a number of these curves drawn and marked with the 
value of f{x, y). If the different values of /(x, y) are chosen 
sufficiently close, so that the curves lie sufficiently close in the 
part of the xy plane that our drawing comprises, we are not only 
able to state the value of /(x, y) at any point on one of the drawn 
curves, but we are also able to interpolate with a certain degree 
of accuracy the value of /(x, y) at a point between two of the 
curves. As a rule it will be convenient to choose equidistant 
values of /(x, y) to facilitate the interpolation of the values 
between. The curves may be regarded as the perpendicular 
projection of certain curves on the surface in space, the inter- 
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sections of the surface by equidistant planes parallel to the 
a:y plane. 

The method is the generalization of the scale-representation 
of a function of one variable. For a relation between i and x 
represented by a curve with i as ordinate and x as abscissa, is 
transformed into a scale representation by perpendicularly 
projecting certain points of the curve onto the axis of x, the 
intersections of the curve by equidistant lines parallel to the axis 
of X and marking them with the value of t. A scale division in 
the case of a function of one variable corresponds to a curve in 
the case of a function of two independent variables. 

This method of representing a function of two independent 
variables by a plane drawing or we might also say of representing 
a surface in space by a plane drawing, is used by naval architects 
to render the form of a ship and by surveyors to render the form 
of the earth's surface and by engineers generally. Let us apply 
the method to a problem of pure mathematics. 

The equation 

Z* + p2 + ? = 

defines z as a function of p and 9. Let us represent this function 

by taking p and q as rectangular coordinates and drawing the 

lines for equidistant values of z. 

For any constant value of z we have a linear equation between 

the variables p and 9, and therefore it is represented by a straight 

line. This line intersects the parallels p = 1 and p = — 1 at 

the points g = — z* — 2 and g = -- ^ + z. Let us calculate 

these values for z = 0; =»= 0,1; =»= 0.2 • • • =*= 1.3 and in this way 

draw the lines corresponding to these values of z as far as they 

lie in a square comprising the values p = — 1 to + 1 and 

g = — 1 to + 1. Fig. 44 shows the result. On this diagram 

we can at once read the roots of any equation of the third degree 

of the form 

z» + pz + g = 0, 

where p and q Ke within the limits — 1 to + 1. For p = 0.4 and 
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9 = — 0.2, for instance, we read z = 0.37, interpolating the value 
of z according to the position of the point hetween the lines 
X = 0,3 and z = 0.4. We also see that there is only one real 
root, for there is only one straight line passing through the point. 



On the left aide of the square there is a triangular-shaped region 
where the straight lines cross each other. To each point within 
this legioQ corresponds an equation with three real roots. For 
example, at the point p= — 0.8 and <? = + 0.2 we read z = 
- 1.00; + 0.28; + 0.72. On the border of this repon two roots 
coincide. 

For values of p and q beyond the linuts — 1 to + I the diagram 
may also be used. We only have to introduce z' = zfrn instead 
of z and to choose m sufficiently loi^e. 

Instead of 

:f + pz+q= 
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we obtain 



or dividing by m*, 



m*z'' + pmsf + 5 = 0, 



or 
where 



»^+ 




= 0, 


«^ + pV + ?'» 


=»o, 


p'- 




9 



By choosing a sufficiently large value of m, p' and 9' can be 
made to lie within the limits — 1 to + 1 so that the roots z' 
may be read on the diagram. Multiplying them by m we 
obtain the roots z of the given equation. 

A function of two independent variables need not be expressed 
in an explicit f orm^ but may be given in the form of an equa- 
tion between three variables 

g{u, V, w) = 0, 

and we may consider any two of them as independent and the 
third as a function of the two. The graphical representation 
may sometimes be greatly facilitated by modifying the method 
described before. The curves for constant values of one of the 
three variables, say w^ are not plotted by taking u and v as 
rectangular coordinates, but they are plotted after introducing 
new variables x and y, x s, function of u and y a function of v and 
making x and y the rectangular coordinates. 

In some cases, for instance, we can succeed by a right choice 
of the functions x = <p(u) and y » }l/(v) in getting straight lines 
for the curves w = const. This will evidently be the case, 
when the equation g{u, v, w) ^ can be brought into the form 

a{w)<p{u) + b(w)4/(v) + c{w) = 0, 

a, b, c being any functions of w, ip any f imction of u and ^ any 
function of 9. 
For introducing 



62 GRAPHICAL METHODS. 

« = ^(v)f y = ^W 

the equation will become 

oa: + 6y + c = 0, 

where a, 6, c are constants for any constant value of w. 

As an example let us consider the relation between the true 
solar time, the height of the sun over the horizon, and the declina- 
tion of the sun for a place of given latitude. Instead of the 
declination of the sun we might also substitute the time of the 
year, as the time of the year is determined by the declination of 
the sun. Our object then is to make a diagram for a place of 
given latitude from which for any time of the year and any 
height of the sun the true solar time may be read. 

In the spherical triangle formed by 
the zenith Z, the north pole P (if we sup- 
pose the place to be on the northern 
hemisphere) and the sun S (Fig. 45), the 
sides are the complements of the decli- 
nation fi, the height h, and the latitude 
ip. The angle t at the pole is the hour 
angle of the sun, which expressed in 
time gives true solar time. 
The equation between these four quantities may be written in 

the form 

sin A = sin ^ sin fi -f- cos tp cos h cos t 

The latitude ^ is to be kept constant, so that f, A, S are the only 
variables. 
Now let us write 

a: = cos f , y ^ Anhf 

so that the equation takes the form 

2^ = sin ^ sin S -f ^ cos ^ cos S. 

When X and y are plotted as rectangular coordinates, we obtain 
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a straight line for any value of S. Let us draw horizontal lines 
for equidistant values of A » to 90^ and vertical lines for equi- 
distant values of < = — 180® to + 180® or expressed in time 
from midnight to midnight (Fig. 46). In order to draw the 




straight lines i = const.^ let us calculate where they intersect 

the vertical lines corresponding to a: = — 1 and x ^ + 1 or 

expressed in time corresponding to midnight and to noon. For 

a; = — 1 we have y = — cos (^ + 8), and for a: = + 1 we have 

y = cos (<p — 8). Let us draw a scale on the vertical a: = — 1 

showing the points y = — cos (<p + S) for equidistant values of 

(^ + 8) and a scale on the vertical a: = + 1, showing the points 

y = cos (^ — J) for equidistant values of ^ — J. The scale is 

the same as the scale for h, with the sole difference that the values 

of ^ — fi are the complements of h and the values of ^ + fi the 

complements of — h. For a latitude of 41®, for instance, we 

have 

For a 

June 21 23.6" 

September 23 and March 21 

December 21 -23.5" 



*» + « 


vS 


64.5'* 


17.6" 


41* 


41* 


17.5" 


64.5* 
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The values of #> + 2 and v— i furnish the intersections with 
the verdcab x = — I and a: = + 1, so that the strai^t lines 
can he drawn corresponding to these days of the year. The two 
outward lines are parallel hut the middle line is steeper. Their 
intersections with the horizontal line k^ show the time of 
sunrise and sunset.' Strictly speaking the straight lines do 
not correspond to certain days. The straight line determined 
by any value of S changes its position continually as S changes 
continually. But the changes of S during one day are scarcdy 
appreciable unless the drawing is on a larger scale. 
If in the equation 

ax + by+e=0 

a and b are independent of to, only e being a function of to, all 
the strught lines to = const, ore parallel. In this case we are 
not obliged to draw the 
strdght lines w = const. 
*• It will suffice to draw a 

line perpendicular to the 
lines v> = const, and a 
scale on it that marks the 
points corresponding to 
equidistant values of to. 
On the drawing we place a 
■♦oo sheet of transparent paper 
^g 47 or celluloid,on which three 

stnught lines are drawn is- 
suing from one point in the direction perpendicular to the u-scale, 
v-scale and w-scale (Fig. 47). If we move the transparent material 
without tumiog it and make the first two lines intersect the u-and-t) 
scale at ^ven points, the w-scale will be intersected at the point 
corresponding to the value of w. This method has the advantage 

' That is to say, the moment when the center of the mm would be seen on 
the horizon, if there were no atmoapheric&l refraction. To t&ke account of 
the refraction, the line A » — 0.6° would have to be conaideied in^JMvl of 
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that we can use the same paper for a great many relations of 
three variables, as we can place a great many scales side by side. 
Or, in the case of one relation only, we may divide the region of 
the values u, v, w into a number of smaller regions and draw three 
scales for each of them, placing all the i^scales or t^-scales or 
ti^scales side by side. The drawing will then have the same 
accuracy as a drawing of very much larger size in which there 
is only one scale for each of the three variables. 

§ 10. Depiction of One Plane on Another Plane. — ^Let us now 
consider two quantities x and y each as a function of two other 
quantities u and v 

y = ^(w, »). 

In order to give a geometrical meaning to this relation between 
two pairs of quantities let us consider x and y as rectangular 
coordinates of a point in a plane and u, n as rectangular coordi- 
nates of a point in another plane. We then have a corre- 
spondence between the two points. When the functions ^(w, i?) 
and ^(t^, v) are defined for the values t«, 9 of a certain region, 
they will furnish for every point w, t of this region a point in 
the Q:y plane. Let us call this a depiction of the ut plane on 
the xy plane. Similarly a function of one variable x == ^(t^) 
might be said to depict the u line on the x line. We may there- 
fore say that the depiction of one plane on another plane is, in 
a certain way, the generalization of the idea of a function of one 
variable. Let us suppose ^(w, u) and ^(w, i) both to have only 
one value for given values of u and t for which they are defined. 
Then there will be only one point in the xy plane corresponding 
to a ^ven point in the wo plane. But to a given point in the 
xy plane there may very well correspond several points in the 
uv plane. 

Let us try to explain this by a graphical representation of the 
depiction of planes on each other. For this purpose we draw 
the curves x = const, and y = const, in the ut plane for equi- 
6 
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distant values of x and y. In the xy plane they correspond to 
equidistant lines parallel to the axis of x and to the axis of y. 
The point of intersection of two lines x= a and y = b corre- 
sponds to the points of intersection of the curves 
p(u, v) = a and iC(% *) — &> 
in the uv plane. If in a certun re^on of the la plane, that 
we consider, they intersect only once there is only one point in 
the region of the uv plane considered and one point in the in/ 
plane corresponding to each other. Fig. 48 shows the depiction 
of part of the us plane on part of the xy plane. We have a net 
of square-shaped meshes in the xy plane and corresponding is a 
net of curvilinear meshes in the uv plane. 

Let us consider the curves x — const, in the ur plane as the 
perpendicular projections of curves of equal hraght on a surface 
extended over that part of the vv plane. From any point P 
of the surface corresponding to the values u, v we proceed an 



infinitely small distance, u changing to u + du, vto v-i- dv and 
xtox + dx, where 



du = cos ads, dv = sin adf, 

where ds ^gnifies the length of the infinitely small line from 
u, to u + du, v+ dv in the ue plane and a the angle its direc- 
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tion forms with the positive axis of x. Let PN be a straight line 
whose projections on the u and v axis are equal to d^/du and 
dipjdv and let us write 

dip ^ dip . 

-=fcosX, -^^fsinX, 

r being the positive length of PN and X the angle between its 
direction and the positive axis of x. Then we have 

d^ . dip 
~ dtl "h Tldv^ rdacos (a — X), 

or 

dx 

T" = r cos (a — X). 

dx/ds measures the steepness of the ascent. It is positive when 
the direction leads upward and negative when it leads downward 
and its value is equal to the tangent of the angle of the ascent. 

From the equation 

dx 

^ = r cos (a — X) 

we see that the ascent is steepest for a ^ \, where dxida ^ r. 
The line PN in the u, v-plane shows the perpendicular projection 
of the line of steepest ascent on the surface x = ip{u, v) and the 
length of PN measured in the same unit of length in which u and 
V are measured is equal to the tangent of the angle of the ascent. 
Let us call the line PN the gradient of the f imction ip{u, v) at the 
point u, V. The direction of the gradient is perpendicular to the 
curve ip(u, v) » const, that passes through the point u, v; for in 
the direction of the curve we have 

dx 
and therefore 

If PN' is the gradient of the function rl/(u, v) at the point u, v, the 
angle between PN and PN' must be equal to the angle formed 
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by the curves x = const, and y = const, that intersect at the 
point u, Vf or equal to its supplement according to the angle of 
intersection that we consider. 

Suppose the gradients PN and PN' do not vanish in any of 
the points in the region of the uv plane that we consider and 
that their length and direction vary as continuous functions of 
u and V. Let us further suppose that the gradient PN' (com- 
ponents: d^jdUy d4//dv) is for the whole region on the left side 
of the gradient PN (components: dip/du, dipfdv), or else for the 
whole region on the right side of the gradient PN, then it fol- 
lows that any one of the curves x = const, and any one of the 
curves y = const, can only intersect once in the region considered. 

This may be shown by considering the directions of the curves 
X = const, and y = const, in the uv plane. Let us consider 
that direction on the curve y = const, in which x increases. If 
this direction deviates from PN the deviation must be less than 
90®, because dx/da and therefore cos (a — X) is positive. Let us 
further consider that direction on the curve x = const, in which 
y increases. If it deviates from the direction of PN' the devia- 
tion must be less than 90®. Let us call these directions the 
direction of x (on the curve y = const.) and the direction of y 
(on the curve x = const.). Now if the gradient PN' is on the 
left of the gradient PN the y direction must also be on the left 
of PN (for if it were on the right of PN being perpendicular to 
PN it would form an obtuse angle with PN') and therefore it 
must be on the left of the x direction (for if it were on the right, 
PN' being perpendicular to the x direction would form an obtuse 
angle with the y direction, which we have seen to be impossible). 
Similarly it may be seen, that if PN' is on the right of PN, the 
direction of y will also be on the right of the direction of x. If 
therefore PN' is on the same side of PN in the whole region 
considered, the direction of y will also be on the same side of the 
direction of x for the whole region considered. This excludes 
the intersection of two curves x = const, and y = const, in more 
than one point. For, suppose there are two points of inter- 
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section and we pass along the curve y ^ const, in the direction of 
X. At the first point of intersection we pass over the curve 
X = const, from the side of smaller values of a; to the side of 
larger values of x. Now if the values of x go on increasing 
as we go along the curve y = const, we evidentiy cannot get 
back to a ciu-ve x ^ const, corresponding to a smaller value of x. 
The only possibility of a second point of intersection would be 
that the direction in which the value of x increases on the curve 
y = const, becomes the opposite^ so that in advancing in the 
same direction in which we came x would decrease again. 

The same holds for the curve 
x^ const. If we pass from one 
point of intersection with a 
curve y — const, along a curve 
X ^ const, to a second point 
of intersection with the same 
curve the only possibility is 
that the direction of y also be- 
comes opposite. This is ex- 
cluded as in contradiction with -pia. 40. 
the direction of y being on the 
same side of the direction of x throughout the whole region (Fig.49) 

It will be useful to look at it from another point of view. Let 
us consider a point A in the uv plane corresponding to the 
values u, t and let us increase u and v by infinitely small positive 
amounts du and dv, so that we get four points ABCD, forming a 
rectangle corresponding to the coordinates. 

A :u,v; B :u + du^v; C :u,v + dv; D :u + du,v+ dt. 

In the xy plane these points are depicted in the points A^ 
By C, Z), the intersections of two curves u and u+ du with two 
curves v and v + dv (Pig. 50). 

The projections of the line AB in the o^ plane on the axes of 
codrdinates are obtained by calculating the changes of x and y 
for a constant value of v and a change du in the value of u 
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dxi^—du, dyi^—du. 

Similarly the projections of ^C are obtained by calculating the » 
changes of x and y for a constant value of u and a change do in the 
value of t 



dv 



dv 



Denoting the lengths of AB and AC by dsi and dst and the angles 
that the directions of AB and AC form with the direction of the 
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positive axis of x (the angles counted in the usual way) by 71 
and 72 we have: 

(2x1 s (bi cos 71, dyi = dsi sin 71 



and 
or 



dxi^ dsi cos 72, dyi^ dsf sin 72, 



d^ 



dsi d^ 



dsi 



7- = cos 7i-T-> -r- = sm 71 -j- 
Bu du du du 



and 



We may call 



6^ 



d»2 dif 



d9% 



--co37»^, -^^ami*^^. 






the scale of depiction at A in the direction AB and 
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av y\dv/ \dv/ 



the scale of depiction at A in the direction AC. It is here under- 
stood that the uv plane is the original, which is depicted on the 
xy plane. If we take it the other way the scales of depiction 
in the directions AB and AC are the reciprocal values du/dsi 
and dv/dsi. 

The area of the parallelogram ABCD in the o^ plane is 

According to the way in which the angles 72 and 71 are defined 
sin (7j — 71) is positive, when the direction AC points to the left 
of the direction AB (assuming the positive axis of y to the left 
of the positive axb of x), and sin (72 — 71) is negative, when AC 
points to the right. Now dudv is equal to the area of the rectangle 
ABCD in the uv plane. Therefore the value of 

dffdyp dip 84/ 
dudv dv du 

is the ratio of the areas ABCD in the two planes and its positive 
or negative sign denotes the relative position of the directions 
AB and AC in the ocy plane. We may call this ratio the scale 
of depiction of areas at the point A. 

dipdtf/ dipdyp 
du dv dv du 

is called the functional determinant of the fimctions <p(u, v) and 
^(w, v). 

We have foimd the scale of depiction of lengths in the direc- 
tions AB and AC. Let us now try to find it in any direction 
whatever. From any point A in the uv plane, whose coordinates 
are u and v, we pass to a point D dose by whose coordinates are 
u + Ati, V + Av. In the ary plane we find the corresponding 
points A and D with coordinates (Fig. 51). 
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^ a: = (p{u, v) jj, x+ Ax'^ ip(u+ Au, v + Av) 
y = \f'(w, ») ' y + Ay = i/'Cu + Au, « + Au) 



We expand according to Taylor's theorem, and writing for 

shortness 

we find 

Ax = ipyjAu + ^,At? + terms of higher order. 

Ay = }pyjAu + \^»A« + terms of higher order. 



I 


\ 




V 




^jrf^' 






Att 






V 




u 


— ^ tt 
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The length of AD and the angle of its direction we denote by 

Ar and a in the uv plane and by As and X in the xy plane. 

The limit of the ratio AsJAr^ to which it tends, when D approaches 

A without changing the direction AD is the scale of depiction 

at the point A in the direction AD. 

Writing 

Au = Ar cos a, 
Av = Ar sin a, 
we obtain 

Ax « (^14 cos a + ipv sin a)Ar + terms of higher order. 
Ay = (\^« cos a + \^, sin a)Ar + terms of higher order. 

Dividing by Ar and letting Ar decrease indefinitely, we have in 

the limit 

dx 

-T- = ^„ cos a + ^ sm a. 
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dy 

-V- = \^i» COS a + \^» sin a. 

For dxjdr and dyldr we may also write dsjdr cos X, dsfdr sin X. 

ds 

'T COS X = ^i» COS a + ^» sin a, 

T"sin X = \^« cos a + \^» sin a. 

These equations show the scale of depiction d^ldr corresponding 
to the different directions X in the Xy y-plane and a in the u^ rh 
plane. 

By introducing complex numbers we can show the connection 
stiU better. 

Let us denote 

2il = ^ + \^»t- 

Multiplying the second of the two equations by i and adding 
both they may be written as one equation in the complex form: 

a = 2i cos a + 2ii sin a. 

The modulus of z is the scale of depiction of the uv plane at the 
point A in the direction a. The angle of z gives the direction in 
the xy plane corresponding to the direction a. For a = we 
have a = ai and for a = 90**, a = 2^. 
Let us substitute 

cos a = rt > sui « ~ or 

and write 

ai + V* » gi — g^A' 
^"^ 2 ' ^^ 2 ' 

so that the expression for a becomes 
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g = a^ -f. te"". 

This suggests a simple geometrical construction of the complex 
numbers z for different values of a. ■ The tenn oe" is represented 
by the points of a cirde described by turning the line that 
represents the complex number a round 
the origin through the angles a = • • • 2x. 
The term fo"^ is represented by the 
points of a circle described by turning 
the line that represents b round the ori- 
gin in the opposite direction through the 
angles a = - ■ - - 2«- (Fig. 52). The 
addition of the two complex numbers 
j^ 52. "^ '^'^^ ^ ** ^^^ *^ value of a b easily 

performed. The points corresponding 
to the complex numbers z describe an ellipse, whose two prind- 
pal axes bisect the angles between a and fa. This is easily seen 
by writing 

a = f,«"<^)', 4 = rjcC***"'. 

ao corresponds to the direction bisecting the angle between a 
and b and oi denotes half the angle between a and 6 (positive or 
negative according to the position of a and b). 

or 

ztf-rf = net—"' + ne -c— »< 

= iXi + ft) cos (a — oi) + (fi — ft) sin (a — ai)t. 

Denoting the coordinates of the complex number ze"*** by ( and ^ 
we have 



and consequently the equation of an ellipse 
El Jt 



(ri+rO*^(ri-ri)«" 
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This ellipse turned round the origin through an angle equal to 
ao gives us the points corresponding to z. The principal axes 
are 2(ri + r^) and 2(ri — ri) (Fig. 53). The construction of 




X 



Fig. 63. 

Fig. 53 is obvious. After plotting Zi and z^ we find z^/i and 
— Zi/i by turning ^422 through a right angle to the right and to 
the left. From these points lines are drawn to Zi. The bisection 
of these lines give a and b. 

The figure shows that in case a and b have the same modulus, 
the triangle — Zt/i, Zi, Z%li becomes equilateral and AZ\ is per- 
pendicular to the line joining — Z^fi and Za/i. In this case AZ\ 
and AZt would have the same or the opposite direction. But as 
Zi = ^u + i^uif ail = ^» + ^vh this would mean that fpu^t "" <pv^u 
= 0. 

The radii of the ellipse (Fig. 53) measured in the unit used 
give the different scales of depiction corresponding to the dif- 
ferent directions in the xy plane. We might also say the ellipse 
is the image in the xy plane of an infinitely small circle in the 
uv plane, magnified in the proportion of the infinitely small radius 
to 1, with its center in A. 

Zx corresponds to a = and Z2 to a = 90® and for a = to 90® 
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Z moves on the ellipse from Z\ to Zs through the shorter way. 
— Z\ corresponds to a == 180® and — Zi to a = 270**. Now we 
have shown above that a positive value of the functional deter- 
minant ipu^p — ipv^u means that Zt is on the positive side of Zi, 
so that in this case Z moves in the positive sense (that is, in the 
direction from the positive axis of x to the positive axis of y) with, 
increasing values of a. With a negative value Z moves in the 
opposite direction. 

Let us now suppose that the curves x = const, and y = const, in 
the uv plane intersect except on a certain ciu^e where their direc- 




tt 
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Xi Xf X^ X^ X( 



tions coincide in the way shown in Fig. 54. On this ciu*ve the 
functional determinant D = <pu^v — <pv^u must vanish because 
the directions of the gradients coincide. Let us see what the 
depiction on the xy plane is like. 

Running along one of the curves y = const., say y = yi, 
toward the curve Z) == we intersect the curves x = X4, Xz, x^ 
until at the point A on the curve a: = jti we reach the curve Z) = 0. 
In the xy plane the corresponding path is a parallel to the axis 
of a; at a distance yi passing through Xi, xz, Xi and reaching a 
point A at Xi. If we now proceed on the curve y = yi in the 
uv plane beyond the curve Z) = 0, we again intersect the curves 
X19 Xz, etc., but in the inverse order. Thus the corresponding 
path in the xy plane does not pass beyond A, but turns back 
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through the same points 0:12, yi; xz, yu etc. The same holds for 
any of the other lines y = const. If we trace the line in the 
acy plane that corresponds to the points in the uv plane, where 
the curves x = const, and y = const, touch, we find the depiction 
of the uv plane only on one side of the curve in the xy plane. 
The other side has no corresponding points u, v. However to 
every point C on this side of the curve, there are two correspond- 
ing points C in the uv plane, one on either side of the curve 
Z) = 0. Imagine two sheets of paper laid on the 0^ plane; let 
them both be cut along the curve AB. Retain only the two 
pieces on this side of the curve and paste them together along 
the curve. The uv plane is in this way depicted on the paper 
in such a way that there is one point and one only on the paper 
corresponding to each point in 
the region of the uv plane con- 
sidered. The curve D = in ^-^, 
the uv plane corresponds to the 
rim where the two pieces of pa- 
per are pasted together. Any 
line straight or curved passing 
over the curve Z) = in the uv 

plane,corresponds to a line running from one of the sheets onto the 
other. It need not change its direction abruptly when it reaches 
the rim and passes onto the other sheet. For it may touch the 
rim in the direction of its tangent. This is actually the rule 
and the abrupt change of direction is the exception. Any line 
LAL (Fig. 55) in the uv plane, whose tangent as it crosses the 
curve D = Bt A does not coincide with the common tangent 
of the curves x = const, and y = const, will correspond to a line 
in the xy plane, that does not change its direction abruptly 
when it touches the rim. 

This is best understood analytically. Let us consider corre- 
sponding directions at the points A in the uv plane and in the 
xy plane. We have seen above that corresponding directions 
(Fig. 56) are connected by the equations 
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ds dx 



COS ^ T" ~ IT ~ ^« ^^s a + ^, sin a, 

^ ds dy 
sinX^ ~ ^ ~ ^" ^^^ a + ^» sin a. 

At the point A we have 

Assuming that the gradients at A do not vanish, so that we 

can write 

ipu^ T cos 7 , ^ = r sin 7, 

^u = r' cos 7', ^, = t' sin 7', 

where r and r ' are positive quantities, the equation ^«^«— ^«^«= 

reduces to sin (7 — 7') = 0, that is, 7 = 7' or 7 = t' + 180**. 

It follows therefore that: 

da 
cos Xt" = r cos (a — 7), 

sin Xt" = / cos (a — 7O = * f' cos (a — 7). 

Consequently for all directions a in the uv plane for which 
cos (a — 7) is not zero, we have 



tgX= =fc-. 
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That is to say, we have in the xy plane only one fixed direction 
X and the opposite corresponding to all the different directions 
OL except only a direction for which cos (a — t) = 0. In the 
latter case, that is, when the direction a is perpendicular to the 
direction 7 of the gradient, t. e., in the direction of the curves 
X = const, and y = const., we have 

cos ^ T" — 0* 

sm X-j-= 0. 

Therefore dsjdf ^ and X remains indeterminate. Any direction 
X for which tg X differs from + r'jr corresponds to a fixed direction 
a = Y + QO^'ora = 7 - 90^ while (fo/dr = 0. 

As the curve Z) = is depicted on the rim of the two sheets 
of paper, all those lines that intersect the curve Z) = in a 
direction different from the direction of the curves x = const, 
and y = const, are depicted in the xy plane as curves having 
their tangent at A in common with the rim. All lines in one of 
the sheets of paper that touch the rim at ^ in a direction differ- 
ent from that of the rim must be the depiction of lines in the uv 
plane that reach A in the direction of the lines x = const, and 
y = const. The scale of depiction is zero in the direction of the 
curves x = const, and y = const. In any other direction a 
we find it different from zero for: 



I =•(,«+/*) C08» («- 7). 

It is a maximum in the direction a = 7 or 7 + 180® perpendicular 
to the curves x = const, and y — const. 

It may help to understand all these details if we discuss an 
example where the depiction of the wo plane on the xy plane 
has a simple geometrical meaning, the planes being ground plan 
and elevation of a curved surface in space. The rim in the 
xy plane is the outline of the surface, the projection of those 
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points where the tangential plane is perpendicular to the plane 

of elevation. 
Suppose a cylinder of circular section cut in two half cylinders 

by a plane through its axis. Suppose one of the half cylinders 

in such a position that its axis 
forms an angle i with the 
ground plan, the plan of ele- 
vation being parallel to its 
axis, Fig. 57. Let us intro- 
duce rectangular coordinates 
u, V in the ground plan and 
rectangular coordinates x, y 
in the plan of elevation. A 
point P on the cylinder is de- 
fined by certain values u, v 
which define its ground plan 
and certain values x, y which 
define its elevation. It is 
easily seen from Fig. 57 that 
we have 




and 



X ^ u 



^ ^ cos 



where a is the radius of the section. Now let us consider the 
elevation of the points P as a depiction of their ground plan. 
The functions ip{u^ t) and ^(u, v) in this case are 

ip{u, f)) = w, 



® cos 



and 



tpu =1, tpv — 0; ^« = tg 8, ^, = — 



cos hVd^ — ^ * 



^u^» — fPti^u = — 



cos 81/^"^^ • 
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Fig. 68. 



The functional determinant vanishes f or v = on the line EF. 
The lines y = const, are the intersections of the cylinder with 
horizontal planes. In the plan of 
elevation they are straight hori- 
zontal lines; in the ground plan 
they are ellipses (Fig. 58). As we 
pass along one of these curves we 
cross the line EF in the ground 
plan but we only touch it in the 
plan of elevation, retracing the hori- 
zontal line back again. The lines 
X = const, are straight lines in both 
planes, but in space they corre- 
spond to ellipses. Again as we 
cross EF in the ground plan we 
only touch it in the plan of eleva- 
tion and retrace the vertical line down again. Any curve on 
the cylinder that crosses EF in a direction not perpendicular to 
the plan of elevation is projected in the plan of elevation with 
EF as its tangent. For the real tangent in space lying in the 
tangential plane of the cylinder can have no other projection, if 
not perpendicular to the plan of elevation. In this latter case 

the projection of the tangent b a point 
and the tangent of the elevation is deter- 
mined by the inclination of the osculatory 
plane. 

There is a particular case to be consid- 
ered, when the ciu^e D = in the uv plane 
coinddes with one of the ciu*ves x ^ const, 
or y = const. (Fig. 59), assuming the gra- 
dients of the functions <p(u, v) and \Ku, v) 
not to vanish at the points of this curve. We have seen that at 
a point where D = the scale of depiction must vanish in the 
directions of the curve x *= const, or y = const. Let the curve 
D = coincide with a line x = const., then it follows that the 




Fxo. 50. 
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length of the depiction of this curve is zero and the depiction 
must be contracted in a point. For the length of the depiction 
of a curve x = const, is given by an integral 



/ 







where dr denotes an element of the ciu^e and ds/dr the scale 

of depiction in the direction of the curve. As ds/dr is zero all 

along the curve the integral must necessarily vanish. 

As an example let us con* 

sider 

X = uv, 

y ^ t. 

The lines x = const, in the uv 
plane are equilateral hyper- 
bolas, the lines y =» const, are 
parallels to the axis of u (Fig. 
60). Along the axis of u we 
have at the same time y = 0, 
ar = and Z)= r = 0. The 
whole axis of ti is depicted in 

the point a; = 0, y = of the xy plane. 
Let us finally consider the case where the scale of depiction 

at any point is the same in all directions, though it need not be 

the same at diiferent points. 
Writing as before 

_d^ dy. _ ds^., 
^^ dr'^ dr^^ dr' 

the connection between the scale of depiction ds/dr and the 
angles X, a determining corresponding directions in the xy plane 
and in the uv plane is given by the equation 

a = Zi cos a + a^ sin a. 
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or 

where 

a = i(2i + ^/i)» b = i(ax — Zi/i). 

In the case where the scale of depiction ds/dr, that is to say, the 

modulus of z, is independent of a, one of the constants a or 6 

must vanish, as we see at once from the construction of z (Fig. 
52). Let us consider the case 6 = 0, 

The complex number a may be written | a \ e^, where | a \ 
denotes the modulus of a and ao the angle. Both may vary 
from point to point, but at every point they have fixed values. 
Consequently we have 

"7" = I a I and X = a + a©. 

That is to say, from an angle a determining a direction in the 
MV plane, we find the angle X determining the corresponding 
direction in the xy plane by the addition of a fixed value ao. 
Any two directions a, a' will therefore form the same angle as 
the corresponding directions X, X' in the xy plane. The same is 
true when o = and z = be"^. The only diiference is that in 
this latter case the direction of z rotates in the opposite sense 
with increasing values of a. 

Analytically depictions of this kind are represented by func- 
tions of complex numbers, 

« + yi = f(u + m) or x+ yi= f(u — w). 

Assuming the function to possess a differential coefficient we have 

dx by. -,, ^ ^ 

dx dy. . 
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and ther^ore either 

2i = aili or si= — %/t. 

HeDce in the first case 

a = iCai + 2i/t) = zi. 6 = ife - ft/O = 

and in the second case 

0=0, b = Zi. 

5 11. Other Method* qf Repretenting Relationi between Three 
Variablea. — The depiction of one plane on another may be used 
to generalize the graphical representation of a function of two 
variables or a relation between three variables, as we prefer 
to say. 

As we have seen before, an equation 

gix, y,z) = 

between three variables x, y, z can he represented by taking x 
and y as rectangular co5rdinates and plotting the curves z = 
const. (Fig. 61) for equidistant val- 
ues of z. Suppose now the xy plane 
to be depicted on another plane. 
The lines x = const.,]/ = const, and 
z = const, will be represented by 
three sets of curves. The fact that 
three values x, y, z satisfy the equa- 
tion g(x, y, z) = Q is shown geo- 
metrically hy the intersection of 
Pio, 61, the three corresponding curves in 

one point. 
Another method for representing cerUun relations between 
three variables u, «, w con^sts in drawing three curves, each 
curve carrying a scale. The values of u, e, w are read each on 
one of the three scales. The relation between three values u, «, 
V} is represented geometrically by the condition that the corre- 
sponding points lie on a straight line (Pig. 62). This method b 
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far more convenient than the one using three sets of curves. It is 
less trouble to place a ruler over two points u, v of two curves 
and read the value w on the scale of the third than to find the 
intersection of two curves u = const, and v = const, among sets 
of others, pick out the curve w = const, that passes through the 
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same point and read the value of w corresponding to it. For we 
must consider that the curves corresponding to certain values 
of u and v are generally not drawn, but must be interpolated and 
so must the curve w = const. It is true that interpolations are 
necessary with both methods, but the interpolation on scales 
like those in Fig. 62 is easily done. 

It must however be understood that while the three sets of 
curves form a perfectly general method for representing any rela- 
tion between three variables, the other method is restricted to cer- 
tain cases. In order to investigate this subject more fully we 
shall have to explain the use of line coordinates. 

When we apply rectangiilar coordinates x, y to define a certain 
jK>int in a plane, we may say that x determines one of a set of 
straight lines (parallel to the axis of ordinates) and y determines 
one of another set of straight lines (parallel to the axis of abscissas) 
and the point is the intersection of the two (Fig. 63, 7). A 
similar method may be used to determine a certain straight line 
in a plane. Let x determine a point on a certain straight line, 
X being its distance from a fixed point A on the line measured 
in a certain unit and counted positive on one side and negative 
on the other. Let y define a point on another straight line 
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parallel to the firsts y being its distance from a fixed point B on 
the line measured in the same way as x. The straight line 
passing through the two points is thus determined by the values 
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X and y and for all possible values of x and y we obtain all the 
straight lines of the plane except those parallel to the lines on 
which X and y are measured. For simplicity we choose AB 
perpendicular to the two lines (Fig. 63, 11). Let us call x and y 
the line coordinates of the line connecting the two points x and 
y in Fig. 63, II ^ in the same way as x and y in Fig. 63, J, are 
called the point coSrdinates of the point where the two lines 
X and y intersect. 

A linear equation between point coordinates 

y = m«+ M 

is the equation of a straight line. That is to say, all the points 
whose coordinates satisfy the equation lie on a certain straight 
line. If, on the other hand, we regard x and y as line coordinates 
we find the analogous theorem: all the straight lines whose 
line coSrdinates satisfy the equation 

y = vfvx + /i 

pass through a certain point. The equation is therefore called 
the equation of the point. 

In order to show this let us first draw the line a; = 0, y = m 
{APO in Fig. 64). If now for any value of x we make AR = x 
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and PQ = mx, the point of intersection of RQ and AP must be 
independent of x, for 

PO mx 

Ark ~ ~ ^» 

AO X 

The ratio PO/AO determines the position of and as it is 
independent of x and the positions of A and P are also inde- 
pendent of X, the same is true for 0. 
For negative values of m, PO and 
AO have opposite directions so that 
lies between ^1 and P. 

For a given point 0, we can find 
the corresponding values of m and /i 
by joining with the points A and 
the point corresponding to x = 1. 
If P and Q are the intersections of 
these lines with the line on which y 
is measured, we have BP = /i and PQ = m. Any point in the 
plane thus leads to an equation 

y^^mx + iJL, 

except the points on the line on which x'ia measured. For 
m = the equation reduces to 

that is, the equation of a point on the line on which y is measured. 

Instead of y = mx + /jl, we might also write x = m'y + ft', 
and go through similar considerations changing the parts of x 
and y. This form does not include the points on the line on 
which y is measured, but it does include the points on the line 
on which x is measured. For these we have m' = 0. 

The general equation of a point in line coordinates is given in 

the form 

ax + by + c= 0, 

from which we may derive either of the first-mentioned forms 
dividing it by a or b. 
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Dividing by c another convenient form is obtained. 



ax . 6g__ 



+ 



— c — c 



or writing 



- + ^ = 1. 

arc »o 

a;o determining the point of intersection of the line BO (Ilg. 64) 
and the online, while ^o determines the point of intersection of 
the line AO with the y-line. 
A curve may be given by an equation 

ai{u)x + bi(u)y + ci{u) = 0, 

in which ai(u), bi(u), ci(u) are functions of a variable u. Any 
value of u furnishes the equation of a certain point and as u 
changes the point describes the curve. Let us suppose the curve 
drawn and a scale marked on it giving the values of u in certain 
intervals sufficiently close to interpolate the values of u be- 
tween them. Two. other curves are in the same way given by 

the equations 

(h(v)x + 6i(r)y + 02(v) = 0, 

ai(w)x + bziw)y + Ci(w) = 0, 

and scales on these ciu^es mark the values of v and w. 

Now we are enabled to formulate the condition which must be 
satisfied by the values u, v, w in order that the three corresponding 
points lie in one straight line. If x and y are the line coordinates 
of the line passing through the three points, x and y must satisfy 
all three equations simultaneously. 

Consequently the determinant of the three equations must 
vanish 

and, vice versa, if the equation between u, v, w may be brought 
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into this form where ai, bi, Ci are any functions of u, 02, Ih, Oi any 

functions of v and <h, bz, c^ any functions of w, we can form the 

equations 

aia: + 6iy + ci = 0, 

(W+lhy + Oi-O, 

a^ + biy + ci^ 0, 

and represent them graphically by curves carrying scales for 
u, V, w. The relation between u, v, w is then equivalent to the 
condition that the corresponding points on the three curves lie 
on a straight line. But it must be remembered that only a 
restricted class of relations can be brought into the required f orm, 
so that the method cannot be applied to any given relation. 
The equation of a point 

ax+hy + c^a 

remains of the same form, when the units of length are changed 
for X and y. If x' denotes the number measiuing the same length 
as the number x but in another unit, the two numbers must have a 
constant ratio equal to the inverse ratio of the two units. There- 
fore, by changing the units independently, we have 

X = Xx', y = My', 

and the equation of the point may be written 

aXx' + hyLy' + c = 0, 
or 

aV + 6y + c = 0, 

where a' ^l<a and V = yh. 

It is sometimes convenient to define the line coordinates in 
another way. Let { and 17 denote rectangular coordinates 
measured in the same unit, then the equation of a straight line 
can be written 

17 = tg ^J + IJo, 

where tp is the angle between the line and the axis of ( and 170, 
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the ordinate of the point of intersection with the axis of 17. 
Now let us call tg ^ and i^o the line coordinates of the straight 
line represented by the equation and let us denote them by x 
and y. Thus the values of x and y define a certain straight line 
and any straight line not parallel to the axis of ordinates may 
be defined in this manner. The condition that a straight line 
Xy y passes through a point £, 17 is expressed by the equation 

or 

y = — fa: + 17. 

If we fix the values of x and y^ all the values i, 17 that satisfy this 
equation represent the points of the straight line x^ y and we 
therefore call it the equation of the straight line. If, on the 
other hand, we fix the values of £ and 17, all the values x, y that 
satisfy the equation represent the straight lines that pass through 
the given point (, 17, and therefore we call it the equation of the 
point. 
The more general form 

or + 6y + c = 
can be reduced to 

a c ' 

It therefore represents the equation of the point, whose rec- 
tangular coordinates are f = a/6 and 17 = — c/6. The case 

where 6 = or 

aa: + c = 

represents the equation of a point infinitely far away in the 
direction ip or the opposite direction ^ + 180®, tp being defined by 

tg^ = a:=-^. 

All the straight lines, whose coordinates x, y satisfy the equation 

aa; + c = 
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coiteapond to the same value of x but to any value of y. That 
is to say, they are all parallel and all the straight lines of this 
directioa belong to them. 

Let us now discuss some of the applications of line coordinates 
to the graphical representation of relations between three 
variables. 

The relation 



may be written in the form 

log « + log » = log w, 
or 

I + y = log w, 
when 

ar = log « and y — log v. 

Let us plot X and y as line co- 

Srdinates on two parallel lines (^g. 

65), with scales for the values ot u no. 6S 

and 0. The equations x = log u 

and y = log r may be regarded as the equations of the points of 

these two scales. The equation 

I + y = log 10 

for any value of to is the equation of a point. It can eanly be 
constructed as the intersection of any lines x, y satisfying its 
equation. For instance, the line x = \og w, y = and the line 
X = 0, y = log w. The first line is found by connecting the 
scale divifflon u = w of the u-scale with the point B, the second 
by connecting the scale divtdon t = w ot the D-scale with the 
point A. If the units of x and y are taken of the same length, the 
point of intersection will lie in the middle between the two lines 
carrying the u and v scales on a line parallel to the two other lines 
and the w-scale will be half the dze of the other two (Fig. 65). 
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log « + log P = log M 

ezpreaaes the condition that the three equations 

a; = log «, y = logv, x+ y = logv 

are satisfied simultaneously by the same values of x and y, that 
is to say, that the three points on the u,v,w scales corresponding 
to the values of u, v, v> lie on the same straight line z, y. 
The more general relation 

« V = w, 
where a and fi are any given values, can be treated in the same 
manner. Thus the pressure and volume of a gas undergoing 
adiabatic changes may be represented. In this case we have 

pp* «= w, 
where p denotes the pressure, v the volume and k and w con- 
stants. 

For a given gas it has a given value, but w depends on the 
quantity of the gas considered. 

We write 

a: = log p, y= log f. 
The rdation then takes the form 
x+ky = \ogfc, 

and represents a point which may be con- 
structed by the intersection of any two 
■Lu straight lines x, y, whose coordinates sat- 
isfy the equation, for instance 

Jl to V 

Pio.e6. a: = logM;, y=0 

and 

1 = 0, 1/ = 7 log w. 
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The first line connects the point B (Fig. 66) with the scale 
division p = w ot the p scale and the second line connects the 
point A with the scale division of the v scale for which y =: A; log w. 
A perpendicular from the point of intersection on AB meets it in 
0' and as the ratio AO'/O'B is 
equal to the ratio of the seg- 
ments on the p and v scales 
log w/k log w = 1/k it is inde- 
pendent of w. All the points 
corresponding to different val- 
ues of w lie on the same par- 
allel to the p and v scales and 
the w scale may be obtained 
by a central projection of the 
p scale on this parallel from 
the center B (Fig. 66). We 
might dispense with the con- 
struction of the w scale as 
long as the straight line for 
the w scale is drawn. For in 
using the diagram we gener- 
ally start with values po, Vo 
and want to find other values 
p, V, for which 
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pVr = Po^O • 
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The straight line connecting the scale divisions p and v intersects 
the w scale at the same point as the straight line connecting the 
scale divisions po and vq, so that we need not know the value of 
P0Oo*. It suffices to mark the point of intersection in order to 
find the value of p, when v is given or the value of v when p is 
given. 
Another example is furnished by the equation 



v^ + xw -}- y = 0. 
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If we regard x and y as line coordinates any value of w determines 
the equation of a point. We plot the curve formed by these 
points with a scale on it indicating the corresponding values of w. 
Any values of x and y determine a straight line whose inter- 
sections with the w scale furnish the roots of the equation. Each 
point of the w scale may be constructed by the intersection of 
two straight lines, whose coordinates x, y satisfy the equation, 
for instance 

a: = 0, y = — v^ and a: = — ir, y = 0.^ 

In Fig. 67 the w scale is shown for the positive values w = to 

w = 2.5. 

In the same manner a diagram for the solution of the cubic 

equation 

ti^ + arw + y = 0. 

or of any equation of the form 

tc^ + ocvf" + y = 
may be constructed. 

§ 12. Relations between Four Variables. — ^The method can be 
generalized for relations between four variables. 
Suppose four variables u, v,w,t€ae connected by the equation 

g{u, V, w, t) = 0, 

and let us assume that for any particular value t ^ to the resulting 

relation between u, v, w can be given by a diagram of the form 

considered consisting of three curves carrying scales for u, v and 

w. Let us further suppose that for other values of t the scales 

for u and v remain the same, but the scale for w changes. Then 

we shall have a set of w scales corresponding to different values 

of t. Connecting the points that correspond to the same value 

of w we obtain a network of curves t = const, and w = const. 

(Fig. 68). Any two values u, v furnish a straight line intersecting 

^ For small values of w, this combination is not good because the angle of 
intersection is small. One might substitute x » 2, y » — to* — 2w for the 
first line. 
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the network of curves. The points ot intersection correspond to 
values of t and w that satisfy the given relation. 
Any relation of the form 

y(u)/((. to) + ^(t)ffft «.) + hit, w) = 

may be represented in this way, ^i>(u) denoting any function of 
u, ilr(e) any function of v and/(^ 
w), g(t, w), k(t, w) any functions 
of t and w. 

In this case we need only in- 
troduce the line coSrdinates z, 
y, writing 

X = v(u), y = Ht). 

We then obtain a linear equation 
between x and y, 

fit, w)x + git, w)y + hit, w) - 0. 

which for any ^ven values of t p,Q. gg. 

and w represents the equation of 

a point. For a given value of t and variable values of w we obtain 
a curve t = const, carrying a scale for w and for a series of values 
of t we obtain a set of curves t = const. Similarly for a given 
value of w and variable values of t the equation furnishes a curve 
w = const, carrying a scale for t and a series of values of v> 
furnishes a set of curves w = const. From any pven values 
of u and V the line coordinates x and y are calculated and the 
points where this straight line defined by x and y intersects 
the network of the curves t = const, and w = const, furnish 
the values t, id that satisfy the relation together with the given 
values of u and v. The relation between the height, azimuth, 
declination of a celestial body and the latitude of the point of 
observation may serve as an example. Let h, a, S denote the 
height, azimuth and declination and v the latitude. The angles 
t/2 — V, H^ — h, t/2 — 3 are the three ades of a spherical 
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triangle PZS (Fig. 69) formed by the pole P, the zenith Z and the 
celestial body S. The azimuth is defined as the supplement of 

z the angle PZS. 

The equation is 

sin 8 = sin ^ sin A — cos ^ cos h cos a. 

We ¥nite 

X = cos a, y = sin 5, 

so that the equation becomes 

y ss sin^sinA — a; cos tp cos A. 

* 

We shall in this case use the second system of line coordinates 
where x is the slope of the line measured by the tangent of the 
angle formed with the axis of gp^ ^o 

abscissas and y is the ordinate 
of the intersection with the axis 
of ordinates. If £, 17 denote the 
rectangular coordinates of the 
pointy the equation of the points 
takes the form 



V'= xi+y or y = 17 — far, 

so that in our case we have 

$ = cos <p cos h, 17 = sin <p sin h. 

The curves tp = const, and h = 
const, can be drawn by means 
of these formulas. It is easily 
seen that they are ellipses and 
that the curves ^ = const, are 
the same as the curves h = const. 
For a definite value of ^ and a 
variable value of h we find 
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COS* ip 



+ 



sin* fp 



= 1, 
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and for a definite value of h and a variable value of ^ 

1^ -!— 2 L ■■'• 



COS* h sin* h 

Any of the ellipses intersects all the others and in this way they 
form a network. A point of intersection of the ellipse ^ = ci 
and the ellipse A = pz also corresponds to the values h = ci and 
^ = pj, as the ellipse ^ = ci is identical with the ellipse A = ci 
and <p ^ Oi identical with h = Ci (Fig. 70). The easiest way to 
find this network consists in drawing the straight lines 

f + 17 = cos (^ - A), 

and perpendicular to them the straight lines 

f- 17 = cos(^+ A), 

for equidistant values of ^ + A and ^ — A. The ellipses run 

diagonally through the rectangular meshes formed by the two 

systems of straight lines. The scales for <p and A are written 

on the axis of coordinates, both scales being available for both 

variables. The scale for d is written on the axis of ordinates 

and is identical with the scale for t and A on this axis. For the 

ordinate corresponding to a given value 5 = c is sin c, and this is 

also the ordinate of the point where the ellipse ^ = c or A = c 

intersects the axis of ordinates. The scale for the azimuth cannot 

be laid down in exactly the same way as that for <p, A and d 

because cos a determines the slope of the straight line x, y. 

Let us draw a parallel to the axis of ordinates through the point 

{ = 1, 17 s= and mark a scale for the azimuth on it, making 

iy = cos a (Fig. 70). A line connecting the origin with any scale 

division of this scale has the slope of the line x = cos a,y ^ sin 5. 

To bring it into the position of the line x, y it must be moved 

parallel to itself, until its point of intersection with the axis of 

ordinates coincides with the scale division b. This suggests 

another way of using the diagram. Let a pencil of rays be 

drawn from the origin to the scale divisions of the azimuth scale 

(Fig. 70), and let it be drawn on a sheet of transparent paper 
8 
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placed over the drawing of the ellipses. For any given valu^ 
of d it is moved up or down as the case may be so that the center 
of the pencil coincides with the scale division B. As long as the 
celestial body does not materially alter its declination the dia- 
gram in this position will enable us to find any of the three 
values <p, h, a from the other two. 

As a second example let us consider the relation between the 
declination 5, the azimuth a, the hour angle < of a celestial body 
and the latitude ip of the point of observation. 

The relation is found by eliminating the height h from the 

equation 

sin 5 = sin ^ sin A — cos fp cos h cos a. 

For this purpose we express sin h and cos h by the other angles 

and substitute these expressions for sin h and cos h. 

We have 

cos h = cos 5 sin f/sin a, 

sin A = sin ^ sin 8 + cos fp cos S cos t 

Substituting these values we find 

sin 3= sin^ ^ sin 5+ sin <p cos tp cos 5 cos ^— cos fp cos 5 sin < ctg a, or 

cos* ^ sin 3 = sin ^ cos <p cos 5 cos t — cos ip cos 5 sint ctg a. 

Dividing by cos* <p cos d we finally obtain 

^ • sin < 

tg = tg ^ cos t ctg a. 

cos tp 

In order to represent this relation graphically we introduce line 

coordinates 

X = ctg a and y = tg 3 

and find 

. . sin< 

y = tg <p cos t z. 

*^ ^ cos ^ 

« 

Let us use the second system of line coordinates. The rec- 
tangular coordinates f, 17 of the point represented by the equation 
are found from it equal to: 
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COS ip 

The curves ip = const, ok ellipses, 



The curves t — const, are hyperbolas, 

ain* i cos* ( 

The ellipses and hyperbolas are confocal, the foci coindding 
with the points { = ± 1, ij = 0, so that the curves intersect at 
right angles. 

The scale for #> may be written on the aids of ordinates at the 
points where it intersects the ellipses. It is identical with the 
scale for h, the ordinate in both cases 
being the tangent of the angle mth the 
only difference that h is negative on 
the negative part of the axis and ip is 
not. The scale for i Toay be written 
on one of the ellipses corresponding to 
the largest value of t» that is to be taken 
account of. This ellipse forms the 
boundary of the diagram, so that 
larger values of *> are not represented. 
Corresponding to the azimuth we draw 
a pencil of rays on a sheet of trans- 
parent paper, which is laid on the draw- 
ing of the curves. Tbe center of the 

/ pendl is placed on the scale division S j.,^ J^ 

f and tbe azimuth is equal to the angles 
that the rays form with the positive direction of the aids of or- 
dinates (Fig. 71). It suffices to draw the curves and the rays 
only on one side of the axis of ordinates. At the apex of the 
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hyperbolas the value of t changes abruptly. The line / = 6* b 
meant to start from the focus £ = 1, 17 == 0. When the center of 
the pencil of rays is in the origin the rays form the asymptotic 
lines of the hyperbolas, a = 16® corresponding tot= 1*, 0=30** 
to < = 2^ and so on. 



CHAPTER III. 

The Graphical Methods of the Differential and 

Integral Calculus. 

§ 13. Graphical Integration. — ^We have shown how the ele- 
mentary mathematical operations of adding, subtracting, multi- 
plying and dividing and the inverse operation of finding the 
root of an equation can be carried out by graphical methods and 
how functions of one or more variables may be represented and 
handled. But the graphical methods would lack generality and 
would be of very limited use, if they were not applicable to the 
infinitesimal operations of differentiation and integration. In- 
deed it is here that they are found of the greatest value. In 
many cases, where the calculus is applied to problems of natural 
science or of engineering, the functions concerned are given in a 
graphical form. Their true analytical structure is not known 
and as a rule an approximation by analytical expressions is not 
easily calculated nor easily handled. In these cases it is of vital 
importance that the operations of the calculus can be performed, 
although the functions are only given graphically. 

Let us begin with integration, because it is easier than differ- 
entiation and of more general application. 

Suppose a function y = f{x) given by a curve whose ordinate is 
y and whose abscissa is z. The problem is to find a curve, whose 
ordinate Y is an integral of the function /(x). 



= I f(x)dx. 



Let us assume the unit of length for the abscissas independent 
of the unit of length for the ordinates. The value of Y measures 
the area between the ordinates corresponding to a and x, the 
curve y == f(x) and the axis of x in imits equal to the rectangle 
formed by the units of x and y. 
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In the simple case where f{x) is a constant the equation 
y ^ /(^) = ^ is represented by a line parallel to the axis of z and 



y ^ \ cdx = c{x — a). 
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Y is the ordinate of a straight line intersecting the axis of z at 
the point x = a. The constant c is the change of Y for an 

increase of x equal to 1. 
If P is the point on the 
axis of a; for a; = — 1 and 
Q the point where the line 
y = c intersects the axis 
of ordinates (Fig. 72) the 
desired line is parallel to 
PQ. It is constructed by 
drawing a parallel to PQ 
through the point x — aon the axis of x (Fig. 72, where a = 0). 
When a given value ci is added, so that the equation becomes 

Y = c (a:— a) + ci 

it amounts to the same as when the straight line is moved in the 
direction of the axis of ordinates through a distance Ci. For 
ar = a we then have F = ci, so that we obtain the line 

y = c(a: — a) + Ci, 

by drawing a parallel to PQ through the point ar = a, y = ci. 

In the second place let us assume that the line y = f{z) consists 
of a number of steps, that is to say, that the function has different 
constant values in a number of intervals x ^ xi to x^, aDi to x^, 
etc., while it changes its value abruptly at 3Di, a^, etc. The 
line presenting the integral 



= ff(x)dx 



does not change its ordinate abruptly. It consists of a con- 
tinuous broken line, whose comers have the abscissas Xt^ xt, etc. 
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The directions of the different parts are found in the way just 
described by the pencil of rays from P to the points a, fi, y, etc. 
(Fig. 73), where the horizontal lines intersect the axis of ordinates. 




Fig. 73. 

To construct the broken line we draw a parallel to Pa through 
the point x ^ Xi (in Fig. 73 Xi is equal to 0) as far as the vertical 
X = Xi. Through the point of intersection with the vertical 
a: = a:^ we draw a parallel to P/3 as far as the vertical x = xz. 
Through the point of intersection with the vertical x = ars we 
draw a parallel to Py and so on. 

Finally let us consider the case of an arbitrary function y = f{x) 
represented by any curve. In order to find the curve 



= £f(x)dx 



we substitute for y = f(x) a function consisting of different 
constant values in different intervals and changing its value 
abruptly when x passes from one interval to the next, so that 
the line representing this function consists of a number of steps 
leading up or down according to the increase or decrease of /(ar). 
These steps are arranged in the following way. The horizontal 
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part A1A2 of the first step (Fig. 73) starts from any point Ai 
of the given curve. The vertical part AtBi and the following 
horizontal part B1B2 are then drawn in such a manner that BiBt 
intersects the curve and that the integral of the given function 
as far as the point of intersection Kb is equal to the integral of the 
stepping line as far as the same point. That is to say, the areas 
between the stepping line and the given curve on both sides of 
the vertical part A2B1 have to be equal. When Kb is fixed the 
right position of AiBi may be found by eye estimate. The eye 
is rather sensitive for differences of small areas. Besides a shift 
of A2B1 to the right or to the left enlarges one area and diminishes 
the other so that even a slight deviation from the correct position 
makes itself felt. In the same way the next step B^CiCt is 
drawn with its vertical part B%Ci in such a position that the 
areas on both sides are equal. The integral of the given curve 
as far as Kc will again have the same value as that of the stepping 
line as far as Kc. And so on for the other steps. The integral 
of the stepping line is constructed in the way shown. It is 
represented by a broken line beginning at the foot of the ordinate 
of Ai, The comers lie on the vertical parts of the steps or 
their prolongations. It is readily seen that the broken line con- 
sists of a series of tangents of the integral curve 



= jr/(a:)(fe; 



and that their points of contact with the integral curve lie on 
the same verticals as the points Ai, Kb, Ke, etc. (In Fig. 73 these 
points are denoted 0, 2, 3, • • • .) That these points lie on the 
integral curve follows from the arrangement of the steps which 
make the integral of the given function at Kb, Kc, * * * equal to the 
integral of the stepping line. Now in the points Ai, Kb, Ke •'• 
the ordinates of the given curve coincide with those of the 
stepping line. Hence both integral lines must for these absdssas 
have the same direction. 

1 In Fig. 73 the lower limit is 0. 
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Having constructed the broken line and marked the points 
2, 3, 4, • • • (Pig. 73), the integral curve is drawn with a curved 
ruler so as to touch the broken line in the points, 0, 2, 3, - - -• 
As the given curve does not change its ordinate abruptly the 
integral curve does not change its direction abruptly. The 
drawing shows how well the integral curve is determined by the 
broken line. There is practically no choice in drawing it any 
other way without violating the conditions. 

The ordinate of the integral curve is measured in the same 
unit as the ordinate of the given curve y = f{x). It may some- 
times be convenient to draw the ordinates of the integral curve 
in a scale different from that of the ordinates of the given curve. 
For instance the value of the integral may become so large that 
measured in the same unit the ordinates of the integral curve 
would pass the boundaries of the drawing board, or else they may 
be so small that their changes cannot be measured with sufficient 
accuracy. In the first case the scale is diminished, in the latter 
case it is enlarged. This is done by altering the position of the 
point P, the center of the pencil of rays that define the directions 
of the broken line. If P approaches the directions Pa, Pp, • • • 
become steeper to the same degree as if keeping P unchanged we 
had increased the ordinates of A1A2, B1B29 - • • in the inverse pro- 
portion of the two distances PO, Hence by diminishing the 
distance PO the ordinates of the resulting broken line are enlarged 
in the inverse proportion. On the other hand, by increasing the 
distance PO the ordinates of the resulting broken line are di- 
minished in the inverse proportion of the distances, because the 
change of the directions Pa, P/3, • • • caused by a longer distance 
PO is the same as if the ordinates of A1A2, B\B%, • • • were di- 
minished in the inverse proportion. The broken line constructed 
by means of the longer distance P'O will therefore be the same as 
if the ordinates of the stepping line were diminished. It therefore 
leads to an integral curve whose ordinates are diminished in the 
same proportion (Fig. 74). 

The graphical integration of 
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is not limited to values x > a. The method is just as well applic- 
able to the continuation of the integral curve for x < a. The 




steps have only to be drawn from right to left. The lower limit 
a determines the point where the integral curve intersects the 
axis of X. 

There is a method for the construction of the vertical parts 
of the steps, which may in some cases be useful, though as a rule 
we may dispense with it and fix their position by estimation. 
Suppose that A and B (Fig. 75) 
are two points where the curve is 
intersected by the horizontal parts 
of two consecutive steps and that 
the curve between A and B is a 
parabola whose axis is parallel to 
the axis of x. The position of the 
Yia. 76. vertical part of the step between A 

andB can be then found bya simple 
construction. Through the center C of the chord AB (Fig. 75) 
draw a parallel CD to the axis of x, D being the point of inter- 
section with the parabola. The vertical part EH of the step in- 
tersects CD in a point whose distance from C is twice the distance 
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from 2>. That this is the right position of EH is shown as soon 
as we can prove that the area ADBGA is equal to the rectangle 
EHBG. The area ADGBA can be divided in two parts, the tri- 
angle ABO and the part ADBCA between the curve and the 
chord. The triangle is equal to the rectangle FIBG, while ADBCA 
is equal to two thirds of the parallelogram MNBA, and hence 
equal to the rectangle EHIF. Both together are therefore equal 
to the rectangle EHBG, and the two areas between the stepping 
line and the curve on both sides of EH are thus equal. 

If the curve between A and B is sup- 
posed to be a parabola with its axis par- 
allel to the axis of ordinates the con- 
struction has to be modified a little. 
Through the center C of the chord AB 
(Fig. 76) draw a vertical line CD as far 
as the parabola. On CD find the point 
K whose distance from C is double the 
distance from D and draw through it a 
parallel to the chord AB. This parallel 
intersects a horizontal line through C at a point £. Then EH 
must pass through i. This may be shown in the following way. 
The area between the parabola ADB and the chord AB is equal 
to two thirds of the parallelogram MNBA, MN being the tan- 
gent to the parabola at the point D. If D' is the point of inter- 
section of NN and the horizontal line through C, we have evi- 
dently 

CL = iCD'. 

Therefore the rectangle EHIF is equal to the area ADB A be- 
tween the parabola and the chord and EHBG is equal to ADGBA, 
Any part of a curve can be approximated by the arc of a 
parabola with suflScient accuracy if the part to be approximated 
is suflBciently small. When the direction of the curve is nowhere 
parallel to the axis of coSrdinates, both kinds of parabolas may 
be used for approximation, those whose axes are parallel to the 
aris of X and those whose axes are parallel to the axis of y. But 
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when the direction in one of the points is horizontal (Fig. 76), 
we can only use those with vertical axes and when the direction 
in one of the points is vertical we can only use those with hori- 
zontal axes. Accordingly we have to use either of the two con- 
structions to find the position of the vertical part of the step. 
Do not draw your steps too small. For, although the difference 
between the broken line and the integral curve becomes smaller, 
the drawing is liable to an accumulation of small errors owing 

to the considerable number 
of corners of the broken 
line and little errors of 
drawing committed at the 
comers. Only practical ex- 
perience enables one to find 
I the size best adapted to 
the method. 

Statical moments of areas 
may be found by a double 
graphical integration. Let us consider the area between the curve 
y = /(^) (Fig. 77), the axis of x and the ordinates corresponding 
to :e = and x = (. The statical moment with respect to the 
vertical through a: =' { is the integral of the products of each 
element ydx and its distance ( — :e from the vertical 
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M^ j^ ({ - x)ydx. 



Let us regard if as a function of ( and differentiate it: 

= 0+ \ ydx. 

That is to say, a graphical integration of the curve y 
beginning at « = furnishes the curve whose ordinate is 

m 

dr 



= f(x) 
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Hence a second integration of this latter curve will furnish 
the curve if as a function of {. As if vanishes f or ( = the 
second integration must also begin at the abscissa x = 0. 




Fio. 78. 

Fig. 78 shows an example. Each ordinate of the curve found 
by the second integration is the statical moment of the area on 
the left side of it with respect to the vertical through this same 
ordinate. The ordinate furthest to the right is the statical 
moment of the whole area with respect to the vertical on the 
right. The statical moment of the whole area with respect to a 
vertical line through any point xi is the integral 



r 



(xi — x)ydx. 



Considered as a function of Xi its diffei;pntial coefficient is 






That is to say, the differential coefficient is independent of xi, 
hence the statical moment is represented by a straight line. As 
its differential coefficient is represented by a horizontal line 
through the last point on the right of the curve 



r 



ydx, 
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the direction of the straight line is found by drawing a line 
through P and through the point of intersection Q of the hori- 
zontal line and the axis of ordinates (Fig. 78). The position of 
the straight line is then determined by the condition that 



X 



(ari — z)ydz 



for OTi == ( is equal to the statical moment 



M{& =J^ {i-x)ydx. 



We have therefore only to draw a parallel to PQ through the 
last point R of the curve for Jf ({) found by the second integration. 
The ordinates of this straight line for any abscissa Xi represent 
the values of 

{zi — z)ydx 



i 



measured in the unit of length of the ordinates. The point of 
intersection E with the axis of z determines the position of the 
vertical in regard to which the statical moment is zero, that is to 
say, the vertical through the center of gravity. 
The moment of inertia of the area 



f 



ydx 



about the axis x =c | is found in a similar way. It is exDiessed 
by the integral 

r= r ({-a:)^<fe. 

Considered as a function of | we find by differentiation 

dT r*d 



= t({ - a:)»yUf + X ^ ({ - «)Vir 



= + 2j ({ - x)yix. 
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That is to say^ the differential coe£Seieiit is equal to double the 
statical moment about the same axis. This holds for every value 
of. {. Hence we obtain |7 as a function of f by integrating 
the curve for -!/({). For { = we have T = 0, so that the curve 
begins on the axis of or at (= 0, 
The integral 

ydx 



J." 



is zero for x ^ a. The curve representing the integral has to 
intersect the axis of 2; at 2; = a (admitting values of 2 > a and 
X < a), and it is there that we begin the construction of the 
broken line. If instead we begin it at the point ar = a, y = c, 
the only difference is that the whole integral ciurve is shifted 
parallel to the axis of ordinates by an amount equal to c upwards 
if c is positive^ downwards if it is negative. But the form of the 
curve remains the same. It is different when this curve is 
integrated a second time. For instead of 



we now integrate 



J ydx 
j ydx + c. 



The ordinate of the integral curve is therefore changed by an 
amount equal to c{x — a) and besides if the second integral curve 
is begun at a: = a, y = ci instead of a: = a, y = the change 

amounts to 

c(x — a) + cu 

so that the difference between the ordinates of the new integral 
curve and the ordinates of the straight line 

y = c(x— a) + ci 

is equal to the ordinates of the first integral curve (Fig. 79). 

This effect of adding a linear function to the ordinates of the 
integral curve is also attained by shifting the pole P upward or 
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downward. For it evidently comes to the same thing whether 
the curve to be integrated is shifted upward by the amount c or 
whether the point P is moved downward by the same amount, so 
that the relative position of P and the curve to be integrated 
is the same as before. Changing the ordinate of P by — c adds 



/ 



/ 



^-"^ --«^r: 



JfdxA-O 
fff(x)dxd» 




ffftx)daodx 
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c(x — a) to the ordinates of the integral curve. c(x — a) is the 
ordinate of a straight line parallel to the straight line from the 
new position of P to the origin. 

By this device of shifting the position of P upward or down- 
ward the integral curve may sometimes be kept within the 
boundaries of the drawing without any reduction of the scale of 
ordinates. A good rule is to choose the ordinate of P about 
equal to the mean ordinate of the curve to be integrated. The 
ordinates of the integral curve will then be nearly the same at 
both ends. The value of the integral 



J ydx 



is equal to the difference between the ordinates of the integral 
curve and the ordinates of a straight line parallel to PO through 
the point of the integral curve whose abscissa is a. 
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When the ordinate of P is accurately equal to the mean ordinate 
of the curve to be integrated for the interval ar = a to 6 the 
ordinates of the integral curve will be accurately the same at the 
two ends. But we do not know the mean ordinate before having 
integrated the curve. 

After having integrated we find the mean ordinate for the 
interval a: = a to 6 by drawing a straight line through P parallel 
to the chord AB of the integral curve, A and B belonging to the 
abscissas x^a and x=b. This 
line intersects the axis of ordi- 
nates at a point whose ordinate 
is the mean ordinate. 

Suppose a beam AB is sup- 
ported at both ends and loaded 
by a load distributed over the 

beam as indicated by Fig. 80. That is to say, the load on dx is 
measured by the area ydx. Let us integrate this curve graph- 
ically, beginning at the point A with P on the line AB. The 
final ordinate at B 
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X 



ydx 



gives the whole load and is therefore equal to the sum of the two 
reactions at A and B that equilibrate the load. Integrating this 
curve again we obtain the curve whose ordinate is equal to 



Y being written for 



J ydx 



The ordinate of this curve at any point x = { represents the 
statical moment of the load between the verticals a; = a and 
X = { about the axis a; = (. Its final ordinate BM, Fig. 81, is 
the moment of the whole load about the point B, and as the reac- 
tions equilibrate the load it must be equal to the moment of the 
9 
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reactions about the same point and therefore opposite to the 
moment of the reaction at A about B. If the reaction at 2I b 
denoted by F« we therefore have 

FaQ> -a)^ C Ydx. 
That is to say. Fa is equal to the mean ordinate of the curve 



= I ydx 



in the interval ar » a to 6. The mean ordinate is found by 
drawing a parallel to AM through P which intersects the vertical 
through A at the point F so that AF ^ Fa. As DB is equal to 
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the sum of the two reactions a horizontal line through F will 
divide BD into the two parts BO = Fa and OD = F». 

Shifting the position of P to P' on the horizontal line FG 
and repeating the integration 

* Ydx, 



£ 



we obtain a curve with equal ordinates at both ends. If we 
begin at A it must end in B. Its ordinates are equal to the 
difference between the ordinates of the chord AM and the curve 
AM (Fig. 81)^ and represent the moment about any point of 
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the beam of all the forces on one side of the point (load and 
reaction). 

The area of a closed curve may be found by integrating over the 
whole boundary. Suppose or = a and 2; == 6 to be the limits of 
the abscissas of the closed curve, the vertical x ^ a touching the 
curve at A and the vertical z = b st B (Fig. 82). By A and B 
the closed curve is cut in two, both parts connecting A and B. 
Let us denote the upper part by y = /i(a:) and the lower part 
by y = /a (a;). The whole area is then equal to the difference 



y 
.'/>' 




^-^ 






l — *-^- 
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or equal to 



£Mx)dx - £ Mx)dx, 

fjl{x)dx + JT f2{x)dx. 



We begin the integral curve over the upper part at the vertical 
X = a at a point E, the ordinate of which is arbitrary, and draw 
the broken line as far as F on the vertical x^ h (Fig. 82). Then 
we integrate back again over the lower part, continuing the 
broken line from F fjo 0. The line EG measured in the unit of 
length set down for the ordinates is equal to the area measured 
in units of area, this unit being a rectangle formed by PO and 
the unit of ordinates. That is to say, the area is equal to the 
area of a rectangle whose sides are PO and EG. 
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The method 13 not limited to the case drawn in Fig. 82, where 
the closed curve intersects any vertical not more than twice. A 
more complicated case is shown in Fig. 83. But in all those cases 




where the object is not to find the integral curve but only to find 
the value of the last ordinate the method, cannot dium to be 
of much use, because it cannot compete with the planimeter. 



For the construction of the broken line we have drawn the 
steps in such a manner that the areas on both sides of the vertical 
part of a step between the curve and the 3tq>pmg line are equal. 
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It would have also been admissible to construct the stepping 
line in such a way that the areas on both sides of the horizontal 
part of a step are equal (Fig. 84). Only the broken line would 
consist of a series of chords instead of a series of tangents of the 
integral curve. The points Ka, Kb, • • • , where the horizontal 
parts of the steps intersect the curve would determine the ab- 
scissas of the points of the integral curve, where its direction is 
parallel to the direction of the broken line. But this forms very 
little help for drawing the integral curve. That is the reason 
why the former method where the broken line consists of a series 
of tangents is to be preferred. However where the object is only 
to find the last ordinate of the integral curve the two methods 
are equivalent. 

§ 14. Graphical DiffereniiaHon, — The graphical differentiation 
of a function represented by a curve is not so satisfactory as the 
graphical integration because 
the values of the differential 
coefficient are generally not 
very well defined by the curve. 
The operation consists in 
drawing tangents to the given 
curve and drawing parallels 
through P to the tangents 
(Fig. 85). The points of in- 
tersection of these parallels 
with the axis of ordinates fur- 
nish the ordinates of the curve representing the derivative. 
The abscissa to each ordinate coincides with the abscissa of the 
pK>int of contact of the corresponding tangent. The principal 
difficulty is to draw the tangent correctly. As a rule it can be 
recommended to draw a tangent of a given direction and then 
mark its point of contact instead of trying to draw the tangent 
for a given point of contact. A method of finding the point of 
contact more accurately than by mere inspection consists in 
drawing a number of chords parallel to the tangent and to 
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bisect them. The points of bisection form a curve that inter- 
sects the given curve at the point of contact (Fig. 86). When a 
number of tangents are drawn, their points of contact marked 
and the points representing the differential coefficient constructed, 

the derivative curve has to be 
drawn through these points. 
This may be done more accur- 
ately by means of the stepping 
line. The horizontal parts of 
the steps pass through the 
points while the vertical parts 
lie in the same vertical as the 
point of intersection of two 
consecutive tangents. The derivative curve connects the iK>ints 
in such a way that the areas between it and the stepping line are 
equal on both sides of the vertical parts of each step. Thus 
the result of the graphical differentiation is exactly the same 
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figure that we get by integration, only the operations are carried 
out in the inverse order. 

A change of the distance T?0 (Fig. 87) changes the ordinates 
of the derivative curve in the same proportion and for the same 
reason that it changes the ordinates of the integral curve when we 
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are integrating, but in the inverse ratio. Any change of the or^ 
dinate of P only shifts the curve up or down by an equal amount, 
so that if we at the same time change the axis of x and draw it 
through the new position of P the ordinates of the curve will 
remain the same and will represent the differential coe£Scient, 

When a function /(ar, y) of two variables is given by a diagram 
showing the curves /(x, y) = const, for equidistant values of 
f{x, y) the partial differential coe£Scients can be found at any 
point xny yo by means of drawing curves whose ordinates represent 
f{x, yo) to the abscissa x oTf{xo, y) to the abscissa y and applying 
the methods explained above. For this purpose a parallel is 
drawn to the axis of x, for instance, through the point a^o, yo 
and at the points where it intersects the curves /(ar, y) = const, 
ordinates are erected representing the values of /(x, ^o) iu any 
convenient scale. A smooth curve is then drawn though the 
points so found and the tangent of the curve at the point Xo 
furnishes the differential coefficient Qfjdx lot a: == aro, y = yo. 

The differential coefficients dfldx, ^I^V are best represented 
graphically by a straight line starting from the point Xy y to 
which the differential coefficients correspond, and of such length 
and direction that its orthogonal projections on the axis of x 
and y are equal to dfjdx and Bfjdy, This line represents the 
gradient of the function /(a;, y) at the point x, y} It is normal 
to the curve /(ar, y) = const, that passes through the point x, y, 
its direction being the direction of steepest ascent. Its length 
measures the slope of the surface z = /(ar, y) in the direction of 
steepest ascent. This is shown by considering the slope in any 
other direction. Let us change x and y by 

r cos a, r sin a 

and consider the corresponding change 

Az = /(a: + ^^ cos a, y + r sin a) — /(ar, y) 

of the function. By Taylor's theorem we can write it 
1 See Chap. II, § 10. 
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— r cos a + T~ f sin a + terms of higher order in r , 
ox ay 

a is the direction from the point x^yto the new point x+ r cos or, 
^ + r sin a and r is the distance of the two points. Dividing 
Az by r and letting r approach to zero we find 

„ Az d/ .df. 

hm — = — cosa+ — sma. 
r dx dy 

This expression measures the slope of the surface z = f{xy) 
in the direction a. Now let us introduce the length I and the 
angle X of the gradient, and write 



Then we have 



f^=/cosX, •f- = isinX. 
dx dy 

-fr COS a + — sin a = / cos (a — X). 
dx dy 



That is to say, the slope in any direction a is proportional to 
cos (a ^ X), it is a maximum in the direction of the gradient 
(a = X) and zero in a direction perpendicular to it and negative 
in all directions that form an obtuse angle with it. When all 
three coordinates are measured in the same unit, the length of 
I measured in this unit is equal to the tangent of the angle of 
steepest ascent. Hence the length of the gradient varies with 
the unit of length. When the unit of length in which the values 
of f{xy) are plotted is kept unaltered, while we change the unit 
of length corresponding to the values x and y^ the length of the 
gradient varies with the square of the unit of length. 

§ 15. Differential Equations of the First Order. — In the problem 
of solving a differential equation of the first order 

I = M ») 

by graphical methods the first question is how to represent 
the differential equation graphically. If x and y are meant to 
be the values of rectangular coordinates, the geometrical meaning 
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of the differential equation is that at every point x, y, where 
f(x, y) is defined^ the equation prescribes a certain direction for 
the curve that satisfies it. Let us suppose curves drawn through 
all those points for which f{x, y) has certain constant values. 
Each curve then corresponds to a certain direction or the opposite 
direction. Let us distinguish the curves by different numbers or 
letters and let us draw a pencil of rays together with the curves 
and mark the rays with the same numbers or letters in such a way 
that each of them shows the direction corresponding to the 
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curve marked with that particular number or letter (Fig. 88). 
Our drawing of course only comprises a certain region in which 
we propose to find the curves satisfying the differential equation. 
It may be that f{xy) is defined beyond the boundaries of our 
drawing. Those regions have to be dealt with separately. 

The graphical representation of the differential equation in 
the region considered consists in the correspondence between 
the curves and the rays. It is important to observe that this 
representation is independent of the system of coordinates by 
means of which we have deduced the curves from the equation 
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We can now introduce any system of coordinates {, 17 and find 
from our drawing the equation 

that is to say^ we can find the value of ^(£, 17) at any point i, ri 
of our drawing. If, for instance, the unit of length is the same 
for ( and 17 we draw a line through the center of the pencil of rays 
in the direction of the positive axis of ( and a line perpendicular to 
it at the distance 1 from the center. The segment on the second 
line between the first line and the point of intersection with one 
of the rays measured in units of length and counted positive in 
the direction of positive 17 furnishes the value of ^({, 17) for all 
the points £, 17 corresponding to that particular ray. In this 
respect the graphical representation of a differential equation 
is superior to the analytical form, in which certain coordinates 
are used and the transformation to another system of coordinates 
requires a certain amount of calculation. 

Now let us try to find the curve through a given point P on 
the curve marked (a) (Fig. 88) that satisfies the differential equa- 
tion. We begin by drawing a series of tangents of a curve 
that is meant to be a first approximation. Through P we draw 
a parallel to the ray (a) as far as the point Q somewhere in the 
middle between the curves (a) and (6). Through Q we draw a 
parallel to the ray (6) as far as 22 somewhere in the middle 
between the curves (6) and (c). Through R we again draw a 
parallel to the ray (c) and so on. The curve touching this 
broken line at the points of intersection with the curves (a), 
(6), • • • is a first approximation. But we need not draw this 
curve. In order to find a better approximation we introduce a 
rectangular system of coordinates x, y, laying the axis of x some- 
what in the mean direction of the broken line. Let us denote 
by yi the function of x that corresponds to the curve forming the 
first approximation. The second approximation y% is then ob- 
tained as an integral curve of /(x, yi), that is, of dyi/dx 
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»^yp+ j fi^fVOdx, 



denoting by Xp, yp, the coordinates of P. For this purpose the 
curve whose ordinates are equal to f(x, yi) or dyijdx has to be con- 
structed first. The values of /(a;, yO are found immediately at 
the points where the first approximation intersects the curve 
(a), (6) • • • by difiFerentiation in the way described above, A 
line is drawn through the center of the pencil of rays parallel to 
the axis of x and a line perpendicular to it at a convenient dis- 
tance from the center. This distance is chosen as the unit of 
length. The points of intersection of this line with the rays de- 
termine segments whose lengths are equal to the values of /(^^ y\) 
on the corresponding curves. These values are plotted as ordi- 
nates to the abscissas of the points where the first approximation 
intersects the curves (a), (6), • • • and a curve 

Y^Kx,y,) 

is drawn (Fig. 88). This curve is integrated graphically begin- 
ning at the point P and the integral curve is a second approxi- 
mation. Again we need not draw the curve. The broken line 
suffices, if we intend to construct a third approximation. In 
this case we have to repeat the foregoing operation. This can 
now be performed much quicker than in the first case because the 
values of /(x, y) on the curves (a), (6), • • • have already been 
constructed and are at our disposal. In order to find the curve 

Y = /(x, y^) 

we have only to shift the same ordinates to new abscissas and 
make these coincide with the abscissas of the points where the 
second approximation intersects the curves (a), (6), •••. The 

curve 

Y = Six, y,) 

is then drawn and integrated graphically, beginning at the point 
P. 
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Suppose now the integral curve did not differ from the second 
approximation, it would mean that 

or that 

that is to say, that y^ satisfies the differential equation. 

If there is a perceptible difference the integral curve represents 
a third approximation. It has been shown by Picard that pro- 
ceeding in this way we find the approximations (under a certain 
condition to be discussed presently) converging to the true solu- 
tion of the differential equation, so that after a certain number 
of operations the error of the approximation must become 
imperceptible. 

Denoting by y» the function of the nth approximation we have 

yn+i = ^p + J /(^' y»)^- 

The true solution with the same initial conditions y = Pp for 
X ^ Xp satisfies the equation 

y^yp+ I fi^f y)dx. 

Hence 

yn+i - y = I [fix, yn) - fix, y)]4x, 

or 

r fi^> yn) - f ix, y) , ., 

y^-i-y^ I - — z—r:, — (»»-»)&. 

•/s^ yn y 

Let us now suppose that the absolute value of 

fix, yn) - fix, y) 
yn-y ' 

for all the values of x, y, yn within the considered region does 



DIFFERENTIAL AND INTEGRAL CALCULUS. 125 

not surpass a certain limit My then it follows that a certain relation 
must exist between the maximum error of yny which we denote by 
en and the maximum error of ^n+i, which we denote by ^n+i. 
The absolute value of the integral not being larger than 

Men \x — Xj, 

( I a: — «» I denoting the absolute value of a? — ar„) we have 

^n+l ^ M\X— Xn\en> 

Hence as long as the distance x^Xn over which the integration 
is performed is so small that 

M\x-Xn\^h<l, 



k being a constant smaller than one, the error of {^n+i cannot be 
larger than a certain fraction of the maximum error of ^n- 
But in the same way it follows that the error of j/n cannot be 
larger than the same fraction of the maximum error of ^n-i* and 

so ouy so that 

e„+i ^ ken ^ i*en-i • • • ^ k^'ei. 

But as ei is a constant and k a constant smaller than one, k^ei 
must be as small as we please for a sufficient large value of n. 
That is to say/ the approximations converge to the true solution. 
M being a given constant the condition of convergence 

Jlfla:-a:p|^jfc<l 

limits the extent of our integration in the direction of the axis of x. 
But it does not limit our progress. From any point P' that we 
have reached with sufficient accuracy we can make a fresh start, 
choosing a new axis of x suited to the new situation. As a 
rule it does not pay to trouble about the value of M and to try 
to find the extent of the convergence by the help of this value. 
The actual construction of the approximations will show clearly 
enough how far to extend the integration. As far as two consecu- 
tive approximations show no difference they represent the true 
curve. 
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Suppose that 

Vn-y 

has the same sign for all values z^ y, yn concerned. Say it is 
negative. Suppose further that y* ^ 2^ is of the same sign for 
the whole extent of the integration 

,^.-,- R'^;':f''' (,.-,w.; 

Jm, yn y 

that is to say, the approximative curve ^n is all on one side of the 
true curve. Then if a; — a:, is positive, ^n^i — y must evidently 
be of the opposite sign from yn — y, or the approximative curve 
y^^i is all on the other side of the true curve from yn. For these 
and all following approximations the true curve must lie between 
two consecutive approximations. If the first approximation yi is 
all on one side of the true curve the theorem holds for any two 
consecutive approximations. This is very convenient for the esti- 
mation of the error. 

In Fig. 88 

f{x, yn) - f(x, y) 

yn-y 

is negative from the point P as far as somewhere near S. The 
first approximation is all on the upper side of the true curve. 
Therefore the second approximation must be below the true 
curve at least as far as somewhere near S. 

When the sign is positive the same theorem holds for negative 
values of X ^ Xp. If the integration has been performed in the 
positive direction of x, it may be a good plan to check the result 
by integrating backwards, starting from a point that has been 
reached and to try if the curve gets back to the first starting 
point. In this direction we profit from the advantage of the 
true curve lying between consecutive approximations and are 
better able to estimate the accuracy of our drawing. 

We have seen that the convergence depends on the maximum 
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absolute value of 

f(Xs Vn) - fix, y) 

yn-y 

for all values of Xy y, pn concerned. In order to find the maximum 
value we may as well consider 

for all values of x, y within the region considered. For if we 
assume df/dy to be a continuous function of y, it follows that 
the quotient of differences 

f(x, yn) - fix, y) 

yn-y 

must be equal to tf/dy taken for the same value of x and a value 
of y between y and yn* This is immediately seen by plotting 
fix, y) as ordinate to the abscissa y tor a fixed value of x. The 
value of the quotient of differences is determined by the slope 
of the chord between the two points of abscissas y and ^n* The 
slope of the chord is equal to the slope of the curve at a certain 
point between the ends of the chord. The value of dfjdy at this 
point is equal to the value of 

fix,yn)-fix,y) 

yn-y 

Now let us consider how the coordinate system may be chosen 
in order to make tf/dy as small as possible and thus obtain the 
best convergence. For this purpose let us investigate how the 
value of df/dy changes at a certain point, when the system of 
coordinates is changed. 

Let us start with a given system of rectangular coordinates i, 
ri with which the differential equation is written 

The direction of the curve satisfying the differential equation 
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forms a certain angle a with the positive axis of £ determined by 

(assuming the coordinates to be measm^ in the same unit). 
Now let us introduce a new system of rectangular coordinates 
X, y connected with the system (, ij by the equations 

a: = { cos w + ij sin w, 
y = — { sin w + 1? cos w, 
which are equivalent to 

5 = a: cos w — y sin w, 
ij = a: sin « + y cos w, 

« being the angle between the positive direction of z and the 
positive direction of ^, counted from £ towards x in the usual way. 
The angle formed by the direction of the curve with the positive 
direction of the axis of a; is a ^ w^ and therefore 

^ = tg (a - ft)) = /(a?, y). 

Consequently we obtain for a given value of ta 

el ^ 1 da 

by " cos^ (a — ft)) dy ' 
or remembering that a is given as a function of £ and ni^ 

^ 

dy cos* 

For simplicity's sake we shall assume that the axis of £ is the 
tangent of the curve ^(£, rf) = const, that passes through the 
given pointy so that da/d^ = 0. 

We then have 

ctf 1 da 

r" ~ — TT \ 7" COS ft), 

dy cos* (a — ft)) di; 

and our object is to find how df/dy varies for different values of 



1 ( da . .da \ 

; ; l — 77Sinft) + — cosft) 1. 

,a - ft)) V d£ drj J 
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u. The value of da/d^ is independent of u; it denotes the value 
of the gradient of a, which we represent by a straight line drawn 
from the origin A (Fig. 89) perpendicular to the curve a = const, 
or p(f , ij) = const. 

It is no restrictioB to assume the value of fla/flij positive; it 
only means that the direction of the positive axis of q is chosen 



in the direction of the gradient. Let us draw the line AB (Fig. 
89) in the direction of the po^tive axis of ( and of the same length 
as the gradient. 

In order to show the values of ^!dy for the different portions 
of the ans of x let us lay off the value of df/dy as an abscissa. 
For instance tot a ~ a, df/dy assumes the value 

da 

— cos a. 

The abscissa corresponding to this value is AB' (fig. 80), the 
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Orthogonal projection of AB on the axis of x. For any other 
portion AC (Fig. 89) correspoDding to some other value of », 
we find da/S^ cos w By orthogonal projection of AB on AC, Tlien 
the division by cos (a — u) furnishes AC and a second division 
by cos (a ~ w) leads to AC. Thus a certain curve can be 
constructed whose polar coordinates are r » ^l&y and u, the 
equation in polar coordinates being 

da umm 

d^ cos* (a — w) 

Id rectangular coSrdlnates (, if the equadon assumes the form 

(cos aj + sin oni)* = ^ C 

This shows that the equation ia a parabola, the axis of which Is 
perpendicular to the direction a. AB' is a chord and the gradient 




A6 is a tangent of the parabola. Bisecting AB' in E^, drawing EK 
perpendicular to AB' as far as the axis of i; and bisecting EK in 
D, we find D the apex of the parabola. The three points A, B', D 
together with the gradient wilt suffice to give us an idea of the 
size and sign of &fli>y for the different positions of the positive 
axis of X. 
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dfl&y vanishes when the axis of x is perpendicular to the curve 
a « const., so that it seems as if this were the most favorable 
position. We must, however, bear in mind that the axis of x 
is kept unaltered for a certain interval of integration. When we 
pass on to other points the axis of x is no longer perpendicular 
to the curve a ^ const, there. The position of the axis of x is 
good when the average value of dffdy is small. In Fig. 90 the 
parabolas are constructed for a number of points on the first 
approximation of a curve satisfying the differential equation. 

If we want to make use of the parabolas to give us the numerical 
values of df/dy the unit of length must also be marked in which 
the coordinates are measured. The numerical value of df/dy 
varies as the unit of length and therefore the length of the line 
representing it must vary as the square of the unit of length. 
But if we draw a line whose length measured in the same unit is 

equal to tttt" * this line would be independent of the unit of 
length. For if Ms the line representing the unit of length and 

Z', V' the lines representing the values cjf/fly and tttt- » df/dy 

oJldy 

would be the ratio V/l and tttt" the ratio V'/l\ hence V = P/r. 

Since V varies as ? with the change of the unit of length V is 
independent of the unit of length. This line V represents the 
limit beyond which the product 

dy 

becomes greater than 1. If df/dy remained the same this would 
mean the limit beyond which the convergence of the process of 

approximation ceases. We might lay off the length of tftt iii 

the different directions in the same way as df/dy has been laid 
off. The result is a curve corresponding, point by point, to the 
parabola, the image of the parabola according to the relation of 
reciprocal radii. But all these preparations as a rule would not 
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pay. It is better to attack the integration at once with an axis 
of X somewhat perpendicular to the curves a = const, as long 
as the direction of the curve forms a considerable angle with 
the curve a = const, and to lose no time in troubling about the 
very best position. The convei^ence will show itself, when the 
operations are carried out. When the angle between the direction 
of the curve that satisfies the differential equation and the curve 
a = const, becomes small the apex of the parabola moves far 
away and when the direction coincides with that of the cur\^e 
a = const, the parabola degenerates into two parallel lines pei^ 
pendicular to the direction of the curve a = const. In this case 
the best position for the axis of x is in the direction of the curve 
a =■ const. Without going into any detailed investigation about 
the best position of the axis of x we can establish the general rule 
not to make the axis of x perpendicular to the direction of the 
curve satisfying the differential equation, that is to say, not to 
make it parallel to the axis of the parabola. But we hardly need 
pronounce this rule. In practice it would enforce its own observ- 
ance, because for that position of the axis of x not only dfldy but 
also /(or, y) are infinite and it would become impossible to plot 
the curve Y = f{x, yi). 

There is another graphical method of integrating a differential 
equation of the first order 

which in some cases may well compete with the first method, 
like the first it is the analogue of a certain numerical method. 
The numerical method starts from given values x, y and cal- 
culates the change of y corresponding to a certain small change 
of X. Let h be the change of x and k the change of ^, so that 
a? + A» y + fc are the coordinates of a point on the curve satisf >dng 
the differential equation and passing through the point x, y. k is 
calculated in the following manner. We calculate in succession 
four values ii, fe, kz, k^ by the following equations — 
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We then form the arithmetica] means 

lc,+ h . h + h 

V = — 2 — « = — 2~' 

and find with a high degree of approximation as long as A ia 
not too large 

i=P+l(«-p)-' 
The new values 

X='x+h, Y=y+k 
are then substituted for x and y and in the same way the coordi- 
nates of a third point are calculated and so on. 

This calculation may be performed graphically in a profitable 
manner, if the function f{x, y) is represented in a way suited to 



Fio. 91. 
'. Kutta, ZeiUehH/t far Malhtmalik und Phytik, Vol. 46, p. 443. 
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the purpose. Let us suppose a number of equidistant parallels 
to the axis of ordinates: x =^ Xo, x ^ Xi, x = Xi, x ^ Xz, • • •. 
Along these lines f(x, y) is a function of y. Let us lay o£P the 
values of /(x, y) as ordinates to the abscissa y, the axis of y being 
taken as the axis of abscissas. We thus obtwi a number of 
curves representing the functions /(xo, y)y fi^^u y)$ f(^f y)i • • •• 
Starting from a point A{xo, yo) on the first vertical x ^ Xo (Fig. 
91) we proceed to a point Bi on the vertical x ^ x^ia the following 
way. By drawing a horizontal line through A we find the 
point A' on the curve representing f(xo, y). Its ordinate is equal 
to f(xo, yo)> Projecting the point A' onto the axis of x we find A'^ 
and draw the line PA'\ P is a point on the negative side of 
the ^-axis and PO is equal to the unit of length by which the 
lines representing f{x, y) are measured. Thus 

OA"IPO = /(xo, ito). 

Now we draw AB\ perpendicular to PA'\ so that if h and h\ 
denote the differences of the coordinates of A and B^ we have 

hilh = OA"IPO, 

h = /(xo, yo)A. 

From Ci the point of intersection of the line ABi and the vertical 
a: = xi we find Ci and Ci" in the same way as we found A' and A" 
from Ay only that Ci is taken in the curve representing the 
values of /(xi, y), and draw the line ABt perpendicular to PCi", 
Denoting the difference of the ordinates of A and Bt by h we have 

or 

fa = /Uo + 2 * yo + 2y*' 

From Ct the point of intersection of the line ABt and the 
vertical x «= xi we find in the same way a point Bt on the vertical 
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x=^ xi and the difference h between the ordinate of ^t and that 
of il is 



Ai =/ f a:o + 2 > yo + 2^ j*. 



From ^8 we pass horizontaUy to Bt on the curve representing 
f(xiy) and vertically down to Bz'\ The line AB^ is then drawn 
perpendicular to PBz", so that the difference k^ between the 
ordinates of B^ and A is 

*4 = f{xo + h,y + h)h. 

The bisection of BtBz and of ^1^4 gives us the points £1 and Et 
and the point B is taken between Ei and E2, so that its distance 
from El is half its distance from E2. The point B is with a high 
degree of approximation a point of the curve that passes through 
A and satisfies the differential equation. 

B is then taken as a new point of departure instead of A, and 
in this manner a series of points of the curve are found. 

In order to get an idea of the accuracy attained the distance 
of the vertical lines is altered. For instance, we may leave out 
the verticals x^ x\ and x » Xz^ and reach the point on the 
vertical x^ x^m one step instead of two. The error of this 
point should then be about sixteen times as large as the error on 
the same vertical reached by two steps, so that the error of the 
latter should be about one-fifteenth of the distance of the two. 
If their distance is not appreciable the smaller steps are evidently 
unnecessarily small. 

The values of /(a:, y) may become so large that an incon- 
veniently small unit of length must be applied to plot them. In 
this case x and y have to change parts and the differential 
equation is written in the form 

dx 1 



dy f(x, y) • 
The values of l/f{x, y) are then plotted for equidistant values of 
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y as ordinates to the abscissa x and the constructions are changed 
accordingly. 

§ 16. Differerdial Equations cf the Second and Higher Orders. — 
Differential equations of the second order may be written in the 
form 

da^ 



^{^^v^ty 



Consequently 



or 



Let us introduce the radius of curvature instead of the second 
differential coefficient. Suppose we pass along a curve that 
satisfies the equation and the direction of our motion is detei^ 
mined by the angle a it forms with the positive axis of x (counted 
in the usual way from the positive axis of x through ninety 
degrees to the po»tive axis of y and so on), s being the length of 
the curve counted from a certain point from which we start. 

We then have 

dy , dx 

^ = tga, ^=cosa. 

<Py _ 1 da _ 1 da 
d^~ co^a dx~ oo^a d$' 

da/ds measures the ^'curvature/' the rate of change of direction 
as we pass along the curve, counted positive when the change 
takes place to the side of greater values of a (if the positive axis 
of X is drawn to the right and the positive axis of y upwards a 
positive value of da/ds means that the path turns to the left). 
Let us count the radius of curvature with the same sign as da/ds 
and let us denote it by p. Then we have 

- = cos* af(x, y, tg a). 

Thus the differential equation of the second order may be said 
to give the radius of curvature as a function of x, y, a, that is to 
say, as a function of place and direction. 




X 
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Let US assume that this function of three variables is repre- 
sented by a diagram^ so that the length and sign of p may quickly 
be obtained for any point and any direction. 

Starting from any given point in any given direction we can 
then approximate the curve satisfying the differential equation 
by a series of circular arcs. Let A (Fig. 92) be the starting point. 
We make MaA perpendicular to the given direction and equal to 
p in length. For positive values of p, Ma must be on the positive 
side of the given direction, for n^ative 
values on the n^ative side. Ma is 
the center of curvature for the curve 
at A. With Ma as center and MaA 
as radius we draw a circular arc AB 
and draw the line BMa> On this line 
or on its production we mark the ,_ 
point Mb at a distance from B equal 
to the value of p that corresponds to 
B and to the direction in which the yiq. 02. 

circular arc reaches B. With Mb as 
center and MbB as radius we draw a circular arc BC and so on« 

The true curve changes its radius of curvature continuously, 
while our approximation changes it abruptly at the points 
A, B, C, "' . The smaller the circular arcs the less will accu- 
rately-drawn circular arcs deviate from the ciu-ve. But it must 
be kept in mind that small errors cannot be avoided, when 
passing from one arc to the next. Hence, if the arcs are taken 
very small so that their number for a given length of curve 
increases unduly, the acciu*acy will not be greater than with 
somewhat longer arcs. The best length cannot well be defined 
mathematically; it must be left to the experience of the draughts- 
man. 

Some advantage may be gained by letting the centers and the 
radii of the circular arcs deviate from the stated values. The 
circular arc AB (Fig. 92) is evidently drawn with too small a 
radius because the radius of the curve increases towards B. If 
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we had taken the radius equal to MhB it would have been too 
large. A better approximation is evidently obtained by making 
the radius of the first circular arc equal to the mean of M^A and 
MhB, and the direction with which it reaches B will also be closer 
to the right direction. 

To facilitate the plotting an instrument may be used consisting 
of a flat ruler with a hole on one end for a pencil or a capillary 
tube or any other device for tracing a line. A straight line 
with a scale is marked along the middle of the ruler and a little 
tripod of sewing needles is placed with one foot on the line and 
two feet on the paper. Thus the pencil traces a circular are. 
When the radius is changed, the ruler is held in its position by 
pressing it against the paper until the tripod is moved to a new 
position. By this device the pencil must continue its path in 
exactly the same direction, while with the use of ordinary com- 
passes it is not easy to avoid a slight break in the ciu-ve at the 
joint of two circular arcs. 

Another method consists in a generalization of the method 
for the graphical solution of a differential equation of the first 
order. 

A differential equation of the second order 



^^{""^y^Tj 



may be written in the form of two simultaneous equations of the 

first order: 

dy 

Si""' 

Let U3 consider the more general form, in which the differential 
coefficients of two functions y, z ot x are given as functions of 
X, y, z: 

£ = f(.x.y,z), 
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dx 



= g{x, y, z). 



We may interpret x,y,z9s the coordinates of a poiDt in space 
and the differential equation as a law establishing a certain 
direction or the opposite at every point in space where /(x, y, z) 
and g{x, y, z) are defined. A curve in space satisfies the dif- 
ferential equation, when it never deviates from the prescribed 
direction. Its projection in the xy plane represents the function 
y and its projection in the xz plane represents the function z. 

Let us represent y and z as ordinates and x as absdssa in the 
same plane with the same sj-stem of coordinates. Any point in 



FiQ. 93. 

space is represented by two points with the same abscissa. The 
functions fix, y, z) and g{x, y, a) we suppose to be given either 
by diagrams or by certain methods of construction or calculation. 
For any point that we have to deal with, the values of /(x, y, z) 
and g{x, y, s) are plotted aa ordinates to the abscissa x, but for 
clearness sake not in the same system of codidinates as y and z, 
but in another system with the same axis of ordinates and an 
axis of X parallel to the first and removed far enough so that the 
drawings in the two systems do not interfere with one another. 
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Starting from a certain point P(xp, yp, Zp) in space we represent 
it by the two points Pi{Xp, yp) and Piixp, Zp) in the first systena 
and the values of f(xp, yp, Zp) and g{Xp, yp, Zp) by the two points 
Ai and At in the second system of coordinates (Pig. 93). The 
points A\ and At determine certain directions MAi, and MAz' 
of the curves x, y and x, z, the point M (Fig. 93) being placed at a 
distance from the axis of ordinates equal to the unit of length by 
which the ordinates representing f(x, y, z) and g(x, y, z) are 
measured. Through Pi and Pt we draw parallels to MAi and 
MAt' as far as Qi and Qt with the coordinates Xq, yq and Xq, Zq. 
With these coordinates the values f(xq, yq, Zq) and g{xq, yq, Zq) 
are determined^ which we represent by the ordinates of the 
points Bu Bt. These points again determine certain directions 
parallel to which the lines QiRi and QtRt are drawn^ etc. In this 
manner we find first approximations yi and zi for the functions 
y and z and corresponding to these approximations we find 
curves representing /(«, yi, Zi) and g(x, yi, zi). These curves are 
now integrated graphically, the integral curve of f(x, y\, zx) 
beginning at Pi and the integral curve of g{x, yi, Zi) at Pt and 
lead to second approximations yt and z^: 



J I f(x, yu sii)dx, 
2f2 = 2p + J g(x, yu zi)dx. 



For these second approximations the values otf{x, y%, z^) and 
g{x, yt, zz) are determined at a number of points along the curves 
X, y% and x, zt suflSciently dose to construct the curves representing 
f{x, yt, Zi) and g(x, y%, z^). By their integration a third approxi- 
mation yz, Zi is obtained 

yz"^ yp+ I f(^f y%y ^)^» 

«8 = 2p + f ?(«> y^y ^)dx, 
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and so on as long as a deviation of an approximation from the 
one before can still be detected. As soon as there is no deviation 
for a certain distance x — Zp the curve represents the true solu- 
tion (as far as the accuracy of the drawing goes). The curve is 
continued by taking its last point as a new starting point for a 
similar operation. 

The distance over which the integral is taken can in general 
not surpass a certain limit where the convergence of the approxi- 
mations ceases. But we are free to make it as small as we please 
and accordingly increase the number of operations to reach a 
given distance. It is evidently not economical to make it too 
small. On the contrary, we shall choose it as large as possible 
without unduly increasing the number of approximations. 

In the case of a diiBPerential equation 



= K^'»'I) 



we have /(or, y, 2) = z, and the curve 2, x is identical with the 
curve representing the values of /(x, y, z). We shall therefore 
draw it only once. 

The proof of the convergence of the approximations is almost 
the same as in the case of the differential equation of the first 
order. 

For the n + 1** approximation we have 

y»H-i - yp-^ \ f(^f Vn, Zn)dx; 2^-1 - Zp+, j g(x, yn, Zn)dx. 

For the true curve that passes through the point Xp, yp, Zp we 
find by integration 



hence 



y = yp+ I f(^f y» 2)^; 2 = Zp + I g(x, y, z)dx; 
Vm-i - y = I [f(p^f »n, 2«) - f{x, y, z)]dx; 

Zn+l - « == J [9(^> yn, Zn) - Q^X, y, z)]dx. 
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Now let US write 

/fa y-i, z„) -fix, y, z) = ,, ,,, (»!. - y) 

yn y 

I /fa y, gn) - /fa y, z) . 

Zn — Z 

and similarly 

. gfa y, gn) - yfa y, z) , 

i ;; ; (g* — g)- 

g|i "" g 

The quotients of differences 

/fa Vn, Zn) - SiXy y, Zn) 

y*»-y 

and the three others are equal to certain values of df/dy, Bffdz, 
dg/dy, dg/dz for values of y, z between y and yn and between z 
and Zn {y, yns z, Zn not excluded). Let us assume that for the 
region of all the values of x, y, z concerned the absolute value of 
^/dy and df/dz, is not greater than Mi, and that of dg/dy and 
dg/dz not greater than M2, and that dn, €n denote the majdmum 
of the absolute values of y — yn and z — z» in the interval 
Xp to X. Then it follows that the absolute values of 

f{Xy Vn, Zn) - /fa y, z) aud g{x, yn, Zn) - jffa y, z) 
are not greater than 

Mi{tn-r^n) and M2(5n+€n). 

Hence for the maximum values of yn+i — y and Zn+i — z, which 
are denoted by dn+i and Cn+i we obtain the limits 

5^-1 ^ -^iC^n + «n) I iC — arp J , €nf 1 ^ Mi^bn +tn)\x—Xj,\^ 

and 

5h-i + €»+i ^ (^1 + Jfi) I a: - arp I (5n + Cn). 

If therefore the interval x — a:, of the integration is so far 

reduced that 

(Ml + M2) I a: - a:p I ^ ]fe< 1, 
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5n4>i + €»4.i is not larger than the fraction k of (5n + cn), but 
from the same reason 

(«» + €n) ^ k(dn^i + e„^i), {dn^i + C,^i) < fc(5„-8 + €„^), CtC.; 

therefore 

That is to say, for a sufficiently large value of n Jn+i and €»+i 
will both become as small as we please. 

As in the case of the differential equation of the first order it is 
not worth while, as a rule, to investigate the convergence for the 
purpose of finding a sufficiently close approximation by graphical 
methods. It is better at once to tackle 'the task of drawing the 
approximations and to repeat the operations until no further 
improvement is obtained. The curve will then satisfy the 
differential equation as far as the graphical methods allow it 
to be recognized. 

When the values of f{x, y, z) or g(x, y, z) become too large 
we can have recourse to the same device that we found useful 
with the differential equation of the first order. Instead of x, 
one of the other two variables y or z may be considered as inde- 
pendent, so that the equations take the form 

^ 1 dz _ g(x, y, z) 

dy "^ f(x, y, z)' dy" f{x, y, z) * 
or 

dx 1 dy _ /(g, y, z) 

dz " g(x, y,z)' dz gix, y, z) ' 

or we may introduce a new system of coordinates x', y', s/ and 
consider the resulting differential equations. 

The second method for the integration of differential equations 
of the first order can also be generalized to include the second 
order. Let us again consider the more general case 

^ = M y» 2), ^ = g(x, y, 2). 
Starting from a point a:, y, 2 the changes of y and z (denoted by 



V Y^.q 2~^ P' = 
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k and J) can be calculated for a small change A of x by the fol- 
lowing formulas analogous to those used for one differential 
equation of the first order: 

h = fix, y, z)h; h = g(z, y, 2)h; 

fe=/(ar+2, y+^f ^+2)*' ^=9\^+2' y^'i* ^+2)*' 

*3=/(a:+2, y+-^> ^+2)*' h=g\x+^, y'^2' ^"^2)*' 
h = fix + h, y+ fe, 2 + k)h; U- g{x+h, y+h, z+ h)h; 

h + *3 ki*+' ki , k + h , _ h+U 

2 ' ^ 2 ' 

and with a high degree of approximation, 

A=P+i(g-p); l^p' + W-pr 

These calculations may be performed graphically. For this 
purpose the functions f{Xy y, z) and g{x, y, z) must be given in 
some handy form. We notice that in oiur formulas the first 
argument assumes the values a;, a: + A/2, a: + A. In the next 
step where a:+A, y+i, z + Z are the coordinates of the starting 
point that play the same part that x, y, z played in the first 
step, we are free to make the change of the first argument the 
same as in the first step, so that in the formulas of the second 
step it assumes the values x + h, x+ ^h, x+ 2h and so on for 
the following steps. All the values of the first argument can 
thus be assumed equidistant. Let us denote these equidistant 
values by 

Xo, Xi, Xi, Xz, • • • . 

The values of /(x, y, z) and g{x, y, z) appear in all our formulas 
only for the constant values 

X = Xo, Xi, Xi, • • • . 

For each of these constants / and g are functions of two inde- 
pendent variables and as such may be represented graphically 
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by drawings giving the curves / = const, and g = const., each 
value of z corresponding to a separate drawing. These drawings 
we must consider as the graphical form in which the diiBPerential 
equations are given. It may of course sometimes be very tire- 
some to translate the analytical form of a differential equation 
into a graphical form, but this trouble ought not to be laid to 
the account of the graphical method. 

The method now is similar to that used for the differential 
equation of the first order, y and z are plotted as ordinates in 
the same system in which x is the abscissa. Equidistant parallels 
to the axis of ordinates are drawn 

X = aro, a? = iti, x ^^ Xi^ etc. 

On the first a: == a^o we mark two points with ordinates yo and ^o, 
and from the drawing that gives the values of fipcuy y, z) and 
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gi^o, y, 2) as functions of y and z we read the values /(aro, yoy zo) 
and g(xo, yof ^o) and draw the lines from xq, yo, and 0:0, zo to the 
points 

a%, yo + h and Xi, Zo + h. 

The intersections of these lines with the parallel x ^ Xi furnishes 

the points 
11 
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xif Wo + "2" and «!, zo + 2 • 
With these ordinates we find from the second drawing the values 

f{xi,yo + '^,zo + ^) and glxuyo+-^,ZQ + ^J, 

and by their help we can draw the lines from xq, yo and Xo, zq 
to the points 

X2,yQ + h and xt, yo + h- 

The intersections of these lines with the line x =^ Xi furnishes the 
points 

Xi, yo+ '2 and xi, 20 + 2 » 
and with these ordinates we find the values 

/ f a?b yo + -g* 2o H- 2 ) ' S' ( «b yo + -g , 2o + 2 j * 
which enable us to draw the lines from Xq, yo and Xq, zq to X2, 

yo "h ^8 and 352, Zq + ^. 

With these two ordinates we find from the third diagram (x = X2) 
the values 

f(xi, yo + Kzo+ k) and g(x2, yo +hfZo+ k), 

which finally enable us to draw the lines from X(^o and 2?oZo to 
a%, yo + h and a^, Zo + /4. 

On the vertical line ;r == x^ we thus obtain four points, ^i, JS2, 
Bz, Ba, corresponding to yo + h, yo + h, yo+ h, yo + h and 
four points, Bi, B2, ft', ft', corresponding to Zo + t, Zo + k, 
Zo + kzo + k (Fig. 94). 

ftft and ftft are bisected by the points Ci and C2; ft 'ft' 
and ft'ft' by the points C/, C2'. Finally CiCt and Ci'Cj' are 
divided into three equal parts and the points B and B' are found 
in the dividing points nearest to Ci and Ci'. 

The same construction is then repeated with B and B' as 
starting points and furnishes two new points on the vertical 
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X =^ Xi and so on. To test the accuracy the construction is 
repeated with intervals of x of double the size. The difference 
in the values of y and of z found for x ^ Xi enables us to estimate 
the errors of the first construction — ^they are about one-fifteenth 
of the observed differences. 

Both methods are without difficulty generalized for the integra- 
tion of differential equations of any order. We can write a 
differential equation of the nth order in the form 



d^x ( dx ^b^\ 



or in the form of n simultaneous equations of the first order 

dx 

dx\ 
dt^^' 






dt 
di 



= f{t, X, Xu X29 '" aTn-l). 



A more general and more symmetrical form is 



dt 
dx\ 



= f%{ti X, Xi, '" Xn^i), 



= fn(t, X, Xi, • • • Xn^i). 



dt 

The functions x, xi, a%, • • • x^^i are then represented as ordinates 
to the abscissa t, so that we have n different ciu^es. When the 
function /(/, x, xi, a^, • • • Xn^i) is given in a handy form, so that 
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its value may be quickly found for any given values of t, x, xi, 
• - - Xn^i, there is no difficulty in constructing n curves whose 
ordinates represent the functions x, xi, x^, • • • Xn^i. Starting 
from given values of t, x, xi, a%, • • • Xf^i we have only to apply 
the same methods that have been explained for the first and the 
second order. 
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